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(57) Abstract 

A computerized method (200) for use in data acquisition and manipulation of two-dimensional patterns in the fields of 
medicine, astronomy, chemistry, biology and biotechnology. The interactive computerized method (200) facilitates matching visu- 
al patterns of polypeptide spots in two-dimensional gel electrophoretograms solubilized into polypeptide constituents that are 
separated by electrophoresis. The computerized method (200) manipulates spot pixel coordinates using staged coordinate trans- 
formation techniques (251) on spot markers and unknown study spots to reduce gel preparation distortions and allows a user to 
produce matching results in an operation (252) that compares the tranformed spot data using either a single reference gel ap- 
proach (255) or a multiple reference gels approach (254) for producing the matching results. The method also includes a spot 
matching verification step (252a) and a step (256) to extract potentially mismatched spots from reported matching results. The us- 
er can also resolve contradictions in resolve step (253) and perform spot matching analysis (270) using isoelectric focusing (PI), 
and molecular weight (MW) dimensional separation data. 
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A COMPUTERIZED METHOD OF MATCHING 
TWO-DIMENSIONAL (E-D) PATTERNS 
Field of the Invention 

This invention relates to computerized methods for data 
5 acquisition and manipulation of two— dimensional patterns in 
the fields of medicine, astronomy, chemistry, biology and 
biotechnology. More particularly, the present invention 
relates to computerized methods for data manipulation of 
visual patterns of polypeptide spots in one— dimension (1— D) 

lO and two-dimensional (E— D) gel electrophoretograms, i.e. 
visual patterns of protein of cells (tissue, or body 
fluids), solubi 1 1 zed into polypeptide constituents that are 
separated by electrophoresis. Even more particularly, the 
present invention relates to interactive computerized 

15 methods of spot data acquisition and pixel coordinate 
manipulation involving coordinate transformation techniques 
and the use of isoelectric focusing (PI), and molecular 
weight (MW) dimensional separation properties during a spot 
matching task. 

EO Description of the Prior Art ; 

In the biotechnology field, two— dimensional 

electrophoresis of polypeptides in polyacry lamide gel is a 
process known for separating the polypeptides, hereinafter 
also referred to as spots, in two dimensions, the first on 

£5 the basis of charge by isoelectric focusing, the second on 
the basis of molecular weight by electrophoresis. The 
resulting two-dimensional gel electrophoretograms from this 
process contain spot patterns that are useful for analyzing 
cell types as well as the genetic metabolic activity of 

30 cells. Since the gels contain thousands of two— dimensional 
spot patterns having geometric characteristics that need to 
be visualized and manipulated for gel comparisons, the need 
for the computational and organizational power of a computer 
system is apparent. The visualization of the spots on the 

35 gel has been made possible by autoradiography or staining 
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techniques, silver-stained electrophoretograms producing 
adequate gel scanning data for analysis, although considered 
less accurate than autoradiograms. The efficiency of the 
results of a spot pattern matching analysis is measured in 
5 percentages of matched spots between gels and depends 
greatly on the accuracy and robustness of the underlying 
computer program that controls data acquisition, spot 
detection and spot pattern data comparison and matching 
tasks. Limitations in the accuracy of data acquisition from 
lO gel images are known to involve several factors including 
film noise, flaws in the gels, which may include localized 
stretching or actual physical breaks in the gels, streaking, 
which are protein complexes formed during gel preparations 
and other positional related distortions caused by added or 
15 subtracted charges that affect mobility of the protein in 
the gel or intensity related distortions caused by faintly 
staining proteins due to heterozygosity characteristics of a 
spot. Certainly, human intervention and introduction of 
artifactual informalities also factors into the equation of 
HO optimizing the matching task. Further, the equipment used 
to create gel images from which the positional and intensity 
data is obtained must also be scrutinized in optimizing the 
spot matching tasks. 

The importance of pursuing and improving quantitative 
analysis of two-dimensional gels has been recognized in the 
biotechnology community and has produced the following prior 
art publications and patents that are of interest in 
considering the present invention. 

For example, an article by P.F. Lemkin et al . entitled 
30 "GELLAB: A Computer System for ED Bel Electrophoresis 
Analysis. II. Pairing Spots", Computers and Biomedical 
Research, Vol. 14, pp. 355-380 <1981>. 

An article by K.P. Vo et al , entitled "Computer 
Analysis of Two-Dimensional Gels", Analytical Biochemistry 
35 Vol. 115, pgs. E58-E71 < 1981 ) . 
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An article by M.M. Skolnick et al . , entitled "Computer 
Programs "for Adapting Two— Dimensional Gels to the Study of* 
Mutation", Clin. Chem. Vol- 28/*, pp. 969-978 (1982). 

An article by M.J. Miller et al, entitled "Computer 
5 Analysis of Two— Dimensional Gels: Semi— Automat ic Matching", 
Clin. Chem. Vol. 26/*, pgs. B67-B75 (1982). 

An article by M.J. Miller et al , entitled "Computer 
Analysis of Two -Dimensional Gels: Automatic Matching", 
Electrophoresis Vol. 5, pgs. 297-303 (1984). 
lO An article by P. Vincens et al , entitled "HERMeSs A 

second generation approach to the automatic analysis o-f two- 
dimensional electrophoresis gels, Part Ills Spot list 
matching", Electrophoresis Vol. 8, pgs. lOO— 107 (1987). 

An article by A.D. Olson et al , entitled "Elsie 4 s 
15 Quantitative Computer Analysis of Sets of Two— Dimensional 
Gel Electrophoretograms" , Analytical Biochemistry Vol. 169, 
pgs. 49-70 (1988). 

M.M. Skolnick et al . , Chapter entitled "An Algorithm 
-for Comparing Two-Dimensional Electrophoretic Gels, with 
SO Particular Reference to the Study of Mutations", in Advances 
in Human Genetics , Vol. 16, Chapter 2, pgs. 55-160, H. 
Harris & K. Hirschom, Editors, (Plenum NY 1986) . 

M.M. Skolnick et al . , in an article entitled "Computer 
Programs for Adapting Two— Dimensional Gels to the Study of 
25 Mutation", Clin. Chem. Vol. 28/*, pp. 969-978 (198E>. 

The prior art patents that relate to the field of the 
present invention concern computerized apparatus and method, 
for spot quantitation, spot analysis, image data 
correspondence, spot detection, image analysis and 
30 comparison, all utilizing pixel x and y coordinate data as 
the primary physical characteristic of the spot or image. 
Included in the above teachings are U.S. Patent Nos. 
4, 638,456, 4,618,937, 4,592,089, 4,590,607, 4,389,670, 
4,812,909, 4,811,218, 4,706,192, 4,741,043 and 4,825,388. 
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Although there have been many computerized advances in 
matching spots in £-D gels, the prior art does not teach 
using staged coordinate trans-formations that -firstly puts 
all of the designated marker spots that highlight a region 
5 or investigative pattern formed by spots in the gels being 
analyzed into a registered relationship with marker spots in 
a designated reference gel, and that secondly, applies the 
resulting marker spot registration relationship data to the 
coordinates o-f the remaining spots bounded by the marker 
lO spots, designated hereinafter as the unknown study spot 
members, to achieve x-y spot localization in the study gels 
having gel preparation distortional effects reduced to 
enable improved accuracy in spot location comparisons and 
improved matching efficiencies. Nor does the prior art 
15 teach the use of PI and MW gel preparation data in spot 
analysis of S-D gels to improve spot matching interpretation 
between two or more E— D gels. 

Therefore, a need is believed to exist for an improved 
interactive computerized method of analyzing two-dimensional 
HO gel electrophoretograms that compensates for gel preparation 
distortional effects using staged coordinate transformation 
techniques prior to comparing and matching polypeptide spots 
contained in the S-D gels. 

Further, a need is seen to exist for an interactive 
25 computerized method for analyzing E-D gels that uses staged 
coordinate transformation techniques and that also utilizes 
PI and MW preparation data that enables a user to combine 
the improved spot localization results associated with the 
transformation technique with the inherent accuracy 
associated with the PI and MW dimensional separation values 
that are at the heart of the electrophoresis process* 

A need is seen to exist for an interactive computerized 
method for analyzing 2— D gels that utilizes the transformed 
spot data resulting from the staged coordinate 
35 transformation technique and that further performs matching 



30 
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verification steps to assure that reported matching spots 
are indeed within target geometrical boundaries. 

A further need is seen to exist for an interactive 
computerized method for analyzing 2-D gels whereby a user 
5 can extract potentially mismatched set of spots to improve 
the accuracy of reported matching results* 

In the broader sense* a need is seen to exist for a 
computerized method for analyzing 2-D patterns in the fields 
of medicine, astronomy, chemistry, biology and biotechnology 

lO that uses staged coordinate transformation, matching 
verification and matching techniques that improve the 
accuracy and efficiency of a 2-D pattern matching task. 
Summary of the Invention 

Accordingly, the primary object of the present 

15 invention is to provide a computerized method of matching 2— 
D image data in the fields of medicine, astronomy, 
chemistry, biology and biotechnology having inherent 
distortional effects resulting from the initial image 
prepar at i on process . 

HO A specific object of the present invention is to 

provide an interactive computerized method of matching 
polypeptide spots in 2— D gels having inherent distortional 
effects resulting from the initial 2-D gel preparation 
process. 

E5 Another specific object of the present invention is to 

provide an interactive computerized method for resolving 
spot matching contradictions in spot clusters that 
inadvertently result during the comparing and matching 
steps. 

30 Yet another specific object of the present invention is 

to provide an interactive computerized method for analyzing 
2-D gels whereby a user can extract potentially mismatched 
set of spots to improve the accuracy of reported matching 
results. 
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A related object of the present invention is to provide 
PI and MW data for the 2— D gel spots being manipulated for 
use during an analysis of output data resulting from the 
comparing and matching subroutines to provide improved 
5 interpretation of the spot matching results* 

According to one aspect of the invention, the foregoing 
objects are accomplished by a computerized method whereby 
each 2— D gel under investigation is scanned to generate an 
initial data file listing each spot's identification, the 
lO gel's name, each spot's original x— y coordinate values » each 
spot's integrated intensity* each spot's area, height and 
width. The user then designates one of the 2-D gels as a 
reference gel and the remaining gels as study gels and 
further designates in each gel, an investigative spot 
15 pattern of interest having a sufficient number of similarly 
positioned spots, generally prominent spots, referred to 
herein as reference m arker spots in the reference gel and as 
study marker spots in the remaining study gels. The 
reference marker spots and the study marker spots form a 
SO boundary for the other spots within the spot patterns which 
are further designated as the unknown reference spots for 
the bounded spots in the reference gel and as unknown study 
spots for the bounded spots in the study gels. The user 
then further includes PI and MW for the designated reference 
25 and study marker spots in each of the gels for later 
utilization. At this point the user has effectively 

generated a modified initial data file for the scanned gels. 

Using the modified initial data files, the user then 
executes several subroutines thats (1) adjust the original 
30 x-y coordinates of the spots in the investigative spot 
patterns in all gels by utilizing a two-stage coordinate 
transformation step and an interpolation step using the 
marker PI and MW data, (2) compare spot coordinates in a 
reference gel with spot coordinates in each study gel 
35 (referred to as a single reference match) to determine 
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potentially matching pairs o-f spots, or that compare spot 
coordinates in one gel with spot coordinates in the other 
gels, each time using one of the gels* including the 
previously designated reference gel, as a reference gel 
5 (referred to as multiple reference match), also to determine 
potentially matching pairs of spots, (3) perform a 
verification step that vectorially manipulates the 
potentially matching spot pair data to acertain whether a 
match spot pair or a non-matching spot pair exists, (*t) 

lO resolve contradicting matching results, (5) compare matched 
spot data results against potentially mismatched spot data 
and potentially correctly matched spot data to improve 
accuracy, efficiency and confidence level of the matching 
results, and <&) that enables a user to manipulate the 

15 various data bases resulting from the operation of the 
subroutines for verifying the results as well as analyzing 
the data from different perspectives. 

Therefore, to the accomplishments of the foregoing 
objects, the invention consists of the foregoing features 

HO hereinafter fully described and particularly pointed out in 
the claims, the accompanying drawings and the following 
disclosure describing in detail the invention, such drawings 
and disclosure illustrating but one of the various ways in 
which the invention may be practiced. 

25 Brief Description of the Drawings ; 

Fig- 1 illustrates a typical hardware configuration for 
performing the interactive computerized method of matching 
protein spot patterns in 2-D gels in accordance with the 
present invention. 

30 Fig. 2 illustrates a flow chart of the operations of a 

software program, that implements the present invention. 

Fig. £a illustrates a flow chart of the adjust 
operation portion of the present invention. 

Figs. 2b, 2c and 2d illustrate a set of flow charts of 

35 the compare and match operation portion of the present 
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invention, that are used to determine potentially matching 
pairs of spots. 

Fig- Se illustrates a flow chart of the extract 
potential mismatch operation portion of the present 
5 invention, that compares the matching results with a file 
that contains matching results that are considered more 
likely to be correct. 

Fig. Sf illustrates a flow chart of the resolve 
contr adit ions operation portion of the present invention 
10 that eliminates ambiguous matching results. 

Fig. 3 illustrates a block diagram overview of a 
computer monitor screen illustrating a menu of the 
interactive computer software program, termed "MATCH WARE" ? 
that implements the present invention. 
15 Fig. ^ is a detailed data flow chart of the operation 

of " M ATCHW ARE " software program that implements the present 
invention. 

Fig. 5 illustrates a pair of gels in an overlay 
relationship to aid in understanding a vector analysis 
80 performed by the second transformation stage of present 
invention. 

Fig. 5a illustrates a planer vector diagram based on 
the graphical gel overlay depicted in Fig. 5. illustrating a 
study spot's nearest marker spot and the next nearest marker 
25 spot. 

Fig. 5b is a diagram illustrating an acceptable 
rectangular area for exact matching spots constructed by 
using the x and y components of vectors depicted in Fig. 5a 
shifted and having a common origin. 
30 Description of the Preferred Embodiment 

Referring first to Fig. 1 where a typical computer 
workstation 100 is illustrated as a means for performing the 
spot matching task of the present invention. System 100 
includes a gel image processing system HO, commercially 
35 available under the trademark "VI BABE" , a personal computer 
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ISO, such as a commercially available "IBM AT", or 
compatible, having a minimum system configuration, shown as 
unit ISOb, o-F 640 kilo-bytes of random access memory, floppy 
disk storage capability* 10 mega-bytes of hard disk storage 
5 and enhanced graphics adapter card. Unit lEOb is coupled to 
a display monitor 120a having color graphic display 
capability, a keyboard 120c, a mouse 120d for cursor 
manipulation. Output from the personal computer ISO can be 
by means of printer 130 for text data/graphical data and a 

10 plotter l^O for graphical display data. In the preferred 
embodiment printer 130 is a laser printer adapted for 
PostScript format, such as a Silentwriter LC 890* 
commercially available from NEC. Also, in the preferred 
embodiment 9 plotter 140 is a plotter compatible with HPGL 

15 graphic plotting format, commercially available as a FACIT 
45SO/4551 . 

Fig- 2 illustrates a flow chart summary of the steps 
involved in performing the 2— D gel matching task in 
accordance with the present invention. In the preferred 

SO embodiment, the steps are performed utilizing a software 
product BOO that will be commercially available under the 
trademark " li ATCHW ARE " from the University of Arizona in 
Tucson, Arizona forthwith filing of this patent application. 
Initially, a gel is scanned to produce gel spot data as 

S5 indicated in steps SIO and SSO. The manner of generating 
the initial spot data information is not the subject of the 
present invention, which generation of data includes the use 
of software programs that are known in the art and are thus 
not describe herein. The data generated includes the spot's 

30 identification, gel name, x-y coordinate values, integrated 
intensity, spot area, spot height and width. The state of 
the art being that the spot listing data available from 
these software programs cannot be used in its raw form to 
produce accurate and efficient spot comparison and spot 

35 matching task due to distortions present in the initial gel 
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preparation. Thus, in accordance with the present 

invention, -further manipulation of the raw data is 
undertaken to compensate -for the distortions. 

The data available at computer ISO is typically a 1024 
5 x lOS^ pixel array of the gel images under investigation. 
Once the spot data is available at computer 1HO, and 
according to the data flow steps 230 and 2*W> shown in Fig. 2 
and further referenced in Figs. 3 and **, the user can 
manipulate the raw data, labeled .1st in Fig. ^ to designate 

lO which gel is to be the reference gel and which gel or gels 
are to be the study gel(s). The user can further manipulate 
the gel image to investigate a windowed spot pattern with 
dominant marker spots and all of the unknown numbered 
reference or study spots, depending upon which gel they are 

15 contained. 

By example, Appendices CI, CH and C3 provide data for 
three gels, labeled "a", "b" and «c", respectively which 
gels are used herein to aid in understanding the present 
invention. Appendix C provides a listing and description of 

20 data files generated by " M ATCHW ARE " during the spot matching 
task, Appendix Cl-<l-6> contains the reference data for the 
gel labeled "a", including hardcopy of the gel's spot image 
and the windowed spot pattern of interest containing 
dominant marker spots labeled m in Appendix Cl-1, and the 

25 corresponding spot data listing in Appendix Cl-(2-6). The 
gel selected as the reference gel is generally determined by 
the operator as representing the best specimen from the 
group of gels, i.e. the one with the least visible 
distortions. Appendices C2-(l-7> and C3-<l-6> contain 

30 information about the two study gels labeled "b" and "c" . 
Included in the Appendices C2 and C3 are hardcopy of the two 
study gel's spot images and the windowed spot patterns of 
interest containing dominant marker spots, also labeled m, 
in Appendices C2-1 and C3-1 and the corresponding spot data 

35 listing C2-<2-7> and C3-(2-6). Collectively, the data is 
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referred to as ElOa -for the reference gel and ElOb and ElOc 
for the two study gels in Fig- ^. The PI and MW data is 
also inputted for later use in the analysis as indicated in 
step 230 in Fig.**. 
5 Once the gels under investigation are set up for 

analysis* step £50, as shown in Fig. 2, is executed to 
selectively perform various subroutines- Included in step 
£50 is a two-stage spot coordinate transformation step E51 
(labeled adjust operation and shown in more detail flow 

lO chart in Fig. Ea and in appendices B(l-23> which comprise a 
source code listing of the adjust subroutine that performs 
the two-stage spot coordinate transformations), a comparison 
and match step E5E (shown in more detail in Figs. 2b, Ec and 
Ed) based on using multiple or single reference gels and 

15 supported by a resolve spot matching contradiction step 253 
( shown in more detail in Fig. Ef > and an extract potential 
mismatched spots step P56 (shown in more detail in Fig. Ee) . 
Appendices A(l — 2) provide a brief description of the 
subroutine steps executed during the adjust operation 

EO illustrated in Fig. Ea. Similarly, Appendices A(3-5) 

provide a brief description of the subroutine steps executed 
during the compare and match operations illustrated in Figs. 
2b, Ec and Ed- Appendix A6 provides a brief description of 
the solve cluster operation step executed during the resolve 

25 contradiction operation illustrated in Fig.Ef. Appendix A7 
provides a brief description of the solve cluster operation 
step executed during the extract potential mismatch 
operation illustrated in Fig.Ee. 

In the adjust operation 251, the two— stage coordinate 

30 transformation includes? (1) performing a first 

transformation step that transforms positional coordinates 
of the reference marker spot members and the unknown 
reference spot members in the reference gel , and each spot 
marker member in the study gel(s) set from the original scan 

35 coordinate system to a new reference coordinate system- 
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This first transformation step results in each member of the 
set of study marker spots being in a registered relationship 
with a corresponding spot marker member of the set of 
reference marker spots. It should be noted that the 

5 coordinates of the spot members in reference gel "a" remain 
the same after the first transformation. <a> performing a 
second transformation step that transforms positional 
coordinates of each of the unknown study spot members from 
the scan coordinate system to the reference coordinate 
lO system, the second transformation step comprises: (a) 
determining an effective range associated with each study 
marker spot, i.e. the distance associated with an influence 
which a dominant marker spot has over the surrounding spots 
which form a spot cluster. The effective range is 
15 calculated for each marker spot by locating the nearest 
marker neighbor to the marker spot and then assigning one- 
half of the distance between the marker and its nearest 
neighbor marker to be its effective range, (b) determining 
an attraction pairing relationship between a particular 
study marker spot member and a particular unknown study spot 
member, this attraction pairing relationship being 
determined by utilizing the effective range of the marker 
spots by a first rule for finding the nearest and second 
nearest markers to the unknown study spot, if this unknown 
25 study spot is within the effective range of the nearest 
marker, then the movement value of this unknown study spot 
into the new reference coordinate system is same as the 
movement value of the nearest marker spot. A second rule 
for determining the pairing relationship can be used if the 
current unknown spot is within the intersection of two times 
the effective range of the nearest and second nearest 
marker, then the movement value of this spot is same as the 
average movement value of the nearest and second nearest 
marker, otherwise the pairing relationship defaults to the 
35 first rule, <c> determining positional coordinates in the 



SO 



30 
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reference coordinate system -For each unknown study spot 
member by adjusting the original scan coordinates for each 
unknown study spot member by a shift amount, or movement 
values, equivalent to the first transformation shift amounts 
5 of the corresponding paired study marker spot member, and 
(d) repeating the pairing and new coordinate location 
determining steps for all unknown study spot members in a 
gel and for all gels being investigated. After the second 
transformation step, the PI and MW values of the unknown 
lO study spot members are determined using interpolation 
techniques based on PI and MW values previously assigned to 
the study marker spot members. Adjusted spot datafiles 851a 
Cx.als), 251b (x.mks) and 251c (x.att) are produced as a 
result of the foregoing two-stage transformation and also 
15 include results of the PI and WW interpolation steps. 
Appendix C4-1 ) contains hardcopy printout of the windowed 
spot pattern and Appendices C**<2-4> contain corresponding 
spot datafiles 251a (a.als), S51b (a.mks) and 251c Ca.att) 
for the reference gel that was manipulated by the 251 adjust 
20 operation, see Appendix C for file descriptors. Appendices 
C5-1 and C6-1 are also hardcopy printout of the windowed 
spot patterns for the two study gels "b" and "c", while 
Appendices C5C2-6) and Co <2-^> contain the corresponding 
spot datafiles 251a Cb.als, c.als), 251b Cb.mks, c.mks) and 
25 251c <b.att, c.att), respectively, resulting from the two- 
stage spot coordinate transformation. 

At this point in the analysis, all the necessary data 
to perform the comparing and matching tasks, denoted in Fig. 
2 as operation 252 and depicted in the flow charts of Figs. 
30 2b, 2c and 2d, is available. The comparing task primarily 
involves taking every identified unknown spot in each of the 
gels and analyzing the recorded data, using as required, the 
spot's old and new x-y coordinate values, integrated 
intensity, spot area, PI and MW data, spot height and width 
35 for comparing against similar data of the other unknown 
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It* 

spots and grouping the results as sets of potential matching 
spots. In the cases involving more than two study gels, the 
user has an option to either conduct a matching exercise 
based on using a single designated reference gel and 
comparing other study gels against and hence following the 
data flow using block 255 shown in Fig. h and in the flow 
chart of Fig. 2b, or conducting a matching exercise 
iteratively designating every gel as a reference gel, i.e. 
multiple reference gels and following the data flow using 
block 25*, also as shown in Fig. * and in the flow chart of 
Fig. Eb. 

Before a matched spots outcome is reported, i.e. based 
on exact matching coordinates, spot size or intensity, the 
present invention also performs a verification step 252a on 
the set of potentially matching spots, regardless of whether 
conducting a matching exercise according to operation step 
25* or 255. The verification inquiry is part of the query 
depicted in the flow chart of Fig. 2d. Verification step 
252a utilizes the potentially matching spot's nearest and 
20 next nearest marker spots, and their original positional 
data, to construct a pair of marker spot vectors for 
juxtaposition comparison with an unknown spot's vector 
formed by joining a first and second one of the potential 
matching unknown spots in two gels, also using the unknown 
25 spot's original positional data. Fig. 5 shows gel Gl and G2 
in an overlay relationship where a first unknown spot G1S1 , 
having nearest marker spot 61M1 and second nearest marker 
spot G1M2, supposedly matching second unknown spot G2S1 , 
having nearest marker spot G2M1 and second nearest marker 
30 spot G2M2. The pair of marker spot vectors comprise a first 
vector VI formed by graphically joining G1M1 to G2M1 and a 
second vector V2 formed by graphically joining G1M2 to G2M2, 
while the unknown spot's vector Ul is formed by graphically 
joining G1S1 to G2S1 , see Fig. 5a. To maintain consistency 
35 in the vector's direction, one of the gels should be 
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designated as containing spots that form the tail of the 
vector and the other as containing spots that "form the head 
of the vectors. The juxtaposition comparison requires 

shifting the three vectors such that their tails are on a 
5 common point to establish whether the head of the unknown 
spot's vector Ul falls within an acceptable rectangular area 
Al formed by minimum and maximum x and y limits, see Fig. 5b 
for construction of the acceptable rectangular area Al , and 
vector Ill's head falling within area Al indicating that 

lO unknown spots G1S1 and GES1 are a matching pair. The size 
of the acceptable rectangular area is based on x and y 
components of the head of vectors VI and VE. The vertical 
boundaries are extended away from the actual x values for 
the VI and VE heads* while the upper y boundary is extended 

15 upward from the y value of the VE vector and the lower 
boundary is extended downward from the y value of the VI 
vector. The x and y amounts extended are user adjustable 
tolerance amounts, typically C5 mi 1 1 i — meters . If the head 
of vector Ul does not fall within the area Al , a no match 

EO situation exists for the particular pair of unknown spots, 
whose vector is being manipulated. 

Assuming that a user wishes to manipulate the gel spot 
data based on multiple reference gel analysis, i.e. 
operation step 25** , then the data flow depicted in Fig. *+ 

E5 would result, see also Fig. 3 for subroutine options 
available. Following the data flow shown in Fig. using 
operation step E5^j the present invention addresses 
ambiguous situations in reported spot matching results (data 
file .mmh shown as E5Eb in Fig. 4) by processing the 

30 affected unknown spots through the resolve contradiction 
subroutine E53, as shown in Fig. and also as shown in the 
flow chart in Fig. Ef. Appendices Dl < 1-3) , DE, D3', DM1- 
3), and D5(l-2) are gel group spot data manipulation results 
in accordance with the present invention wherein: Appendix 

35 Dl is a listing of datafile 25Eb (abc.mmh) based on multiple 
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reference gels, i.e. results based on each gel being used as 
a reference in the comparison task and includes matched and 
unmatched spots be-fore resolving any contradictions 
contained in the data -for the three gels under 
investigation. Appendix D2 is a listing of datafile 250c of 
a spot cluster exhibiting contradicting matching results 
which were extracted from datafile 252b and resolved by 
resolve contradiction subroutine 253, i.e. a spot(s) that is 
(are) found matching other spots that logically cannot be 
explained f for example one spot in one gel matched two spots 
and these two spots are in the same gel, (datafile 250c is 
an empty file upon resolving all contradictory clusters). 
Appendix D3 is a listing of datafile 250d of unmatched spots 
(termed "unique") as found in the three gels under 
15 investigation. Appendix D<*(l-3) is a listing of datafile 
250b of exact matched spots as determined by the present 
invention. Appendix D5( 1-2) is a composite listing of 

datafile 250a of a pseud o gel, each member of the composite 
listing representing a matrix of matching spots having a 
20 plurality of rows, each row being identical to the 
corresponding member in the composite listing. In the 250a 
listing, negative record label represents matching spots, 
(two matching spots or three matching spots), while positive 
record label represents unmatched spots, as found in the 
25 three gels under investigation. 

Referring now to Figs. 3 and 4, Appendix El, is a 
listing of datafile 252c (abc-smh) similar to datafile 252b 
except generated using a designated single gel as a 
reference. Appendix El is generated by operation step 255 
30 that is based on using a designated single gel as a 
reference gel and comparing it against each of the other 
study gels before generating the results- Datafile 252c has 
the characteristic that matched and unmatched spots are 
reported without any matching contradictions which normally 
35 results when the matching is based on a multiplicity of 
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reference gels. It has been observed that some o-f the 
matching spots results include mis— matched spot data and has 
led to reporting more matched spots in matching results 
listings* such as in Appendix El, than actually exist. 
5 Faced with this potential error in matching results using a 
single gel as a reference gel* the user has the option of 
further manipulating the data using step 256 to extract 
potential mismatching spots to improve the accuracy of the 
match listings* see Fig. and also the flow chart depicted 

lO in Fig. He. Step 256 allows the user to compare the 
datafile 252c against a composite datafile 250a, see 
Appendix D5(l-2), abc.cmp, (generated in background during 
operation of step 255 and contains the same data as if 
generated by operation 25**, 253)* which datafile abc.cmp 

15 contains matching results that are considered more likely to 
be correct* Comparison step 256 then generates potentially 
mis-matched spots datafile 256a* see Appendix E2, abc.pmm, 
and a potentially correctly matched spot datafile 256b, see 
Appendix E3* abc.pcm. In generating datafile 256a and 256b * 

20 it is noted that in comparing datafiles 252c against 250a 
that any potentially correct matching set of spots would be 
found in both datafile 250a and in 252c* while any potential 
mis— matched set of spots would be found only in datafile 
252c* but not in 250a. Further, after step 256* it is noted 

25 that any potentially correct matching set of spots would be 
mostly found in datafile 256b and least in 256a* while any 
potential mis— matched set of spots would be found mostly in 
datafile 256a and least in 256b. Composite datafile 250a is 
given greater credibility because the matched set of spots 

30 found in datafile 250a are determined firstly as a matrix of 
matching spots using each gel as a reference gel and then 
repeated according to the number of gels under study = 
Having processed spot data through operation step 255* the 
user can then utilize datafile 256a containing the list of 

35 potentially mis-matched sets of spots as a means of 
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determining whether- a particular set of matching spots being 
examined are part of the potentially mis— matched group of 
spot and thereby gain a higher confidence level about 
reported matching results, i.e. the matching results found 
in datafile 256b or 250b . The spot matching result listed 
in datafile 256b have an improved accuracy over those listed 
in datafile E5Sc in that the potential mis-matching spots 
have been sorted into datafile 256a. 

Referring back to Fig. 2, the results from operation 
250 may be reviewed in a variety of plotting options 260, 
see also Figs. 3 and and may be repeatedly performed 

after analysis step 270 until a user has thoroughly 
understood the matching results, which step 270 takes 
advantage of having the PI and MW data to reinforce matching 
results based on comparisons in the new x and y coordinate 
system, and may be repeated using the data base information 
and may be based on inquiry 280, 281 that considers whether 
any further contradiction 253a exists, or whether new 
markers need to be considered, such as by looping back 
through marker designation step 230 and bypassing step B^O. 

Therefore, while the present invention has been shown 
and described herein in what is believed to be the most 
practical and preferred embodiments, it is recognized that 
departures can be made therefore within the scope of the 
invention, which is therefore not to be limited to the 
details disclosed therein but is to be accorded the full 
scope of the claims so as to embrace any and all equivalent 
methods. 
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ADJUST OPERATION 



ReadMatchCfg 



ReadGetCfg 

ReadListFiles 
. AssignMarkers 



SetMarkersMove 



SetRe fNewCoor 



AdjustXY 



Adjust 



Get the gel names (i.e. file names) from 
match. cfg file. These files will be used to 
compare each other later. The first gel name 
will be the reference gel, the others are 
study gels. 

Get the PI and MW values of each marker. The 
PI vand MW values stored in get. cfg are input 
by user via marker_entry program. 

Get all the spot information of each gel. 

Aligns all the markers of study gels to the 
markers of the reference gel . After 
alignment , all markers have the same X-Y 
value. The X-Y values of the reference gel 
are not changed. Calculates the effective 
range of each marker. 

Find out the movement value of each marker in 
the study gels. The movement value is got 
from the difference of old position and 
aligned position. 

Set the new coordinate of the reference gel. 
(This new coordinate is same as old 
coordinate) 

Kicks out all the spots outside the window 
formed by markers . Call Adjust subroutine to 
adjust all the spots in each gel. 

Adjusts every unknown spot in the study gel 
by the following rules: 

1 . Find out the nearest and second nearest 
markers to this unknown spot. If this 
unknown spot is within the effective 
range of the nearest marker, then the 
movement value of this spot is same as 
the movement value of the nearest mark. 

2. If rule 1 is not satisfied, use this 
rule. If the current unknown spot is 
within the intersection of two (2) times 
effective range of the nearest and 
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second nearest marker, then the movement 
value of this spot is same as the 
average movement value of the nearest 
and second nearest marker. 
If rule 1 and rule 2 are not satisfied, 
then use rule 1 to get the movement 
value of the unknown spot. 



Assigns PI and MW values of the markers in 
each marker list . 



CalMkBoundarylnfo 



GetPiMw 



PrintAttMks 



Calculates the left, right, top, and bottom 
boundary delta x, y, pi, and MW for the using 
in GetPiMw subroutine. 

Calculates the PI and MW values of each spot 
in each gel by using linear interpolation. 
The PI and MW values of each marker are input 
by user. 

Writes out the adjusted information and PI, 
MW value to the files. These files are 
" .als", ".att", and ".mksV The " .att" and 
".mks" will be used in the comparison step. 



NOTE; The algorithm to get the 

marker : 

Find out its nearest marker, 
between this marker and its 
effective range. 



effective range of each 

Get half of the distance 
nearest marker to be its 
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COMPARISON OPERATION 



PCT/US91/03620 



OpenMatchCfg Gets the gel names (i.e. file names) from 

match. cfg file. These files will be used to 
compare each other later. The first gel name 
will be the reference gel, the others are 
study gels. 

ReadMksAtt Gets the marker information of each gel in 

.mks and the spot information of each gel in 
. att . 



Initialize 
Match26els* 

Printout 



Gives initial values to the internal 
parameters of the comparison routine. 

There is an iterative loop to call this 
routine to do the comparison in any pair of 
gels . 

In case that the match type is single 
reference match, writes the comparative 
result to the file (.smh), otherwise, writes 
the result to the file (.mmh). The ".smh" 
will be used in the extract mismatch; the 
".mmh" will be used in the resolve cluster 
step. 



NOTE: The asterisk (*) note means that the routine will be 

described detailedly later. 
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Matcli2Gels description 



Finds out the shift of markers between two 
gels. 



InitMatch 



Gives initial value of tolerance, candidate 
range, minimum neighbor match number, and 
area threshold. 



ClearMatch 



Finds out all the spot pair which are very 
close. All of these spot pairs are clear 
matched. Takes these spot pairs from the 
unmatched group. 



Keep a fixed times loop to the following routines: 



CalNewParameter 



Calculates the new parameters for criteria of 
comparison. If the spot pair satisfies these 
criteria, we say it is matched. 



For eac h spot in first gel, do the following routines 
PindCandidate 

Mat ch2 Spot* 
ManyToOne 



Finds out the match candidates (the spots in 
second gel which are fallen within the 
candidate range) for the current spot in the 
first gel. 

Compares pairs which are fallen within the 
candidate range . 

Solves the status that two spots in the first 
gel matches with one spot in the second gel 
using some creteria. The scores of the 
previous and current matched pair depend on 
area ratio and shift length. 
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SecondChance 
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Mat ch2 Spot description 

Finds out the number of neighbors of the 
first gel which are matches with the 
neighbors of second gel. If the number is 
greater than or equal to the minimum neighbor 
match number, then call I sCorrectMove routine 
to double check matching, otherwise, call 
SecondChance routine to give another chance 
to this spot pair. The reason of the 
existing of IsCorrectMove and SecondChance is 
used to make up the defect of MatchPattern 
routine. In the case that the spot pair 
match is confirmed by IsCorrectMove, they are 
matched and taken from unmatched group, 
otherwise they are not matched. Same case 
will happen when the program calls 
SecondChance routine . 

Solves the status that one spot in the first 
gel matches with two spots in the second gel 
using some creteria. The scores of the 
previous and current matched pair depend on 
area ratio and shift length. 

Finds out whether the spot shift is similar 
to the shift of it's near marker (includes 
shift length and shift angle) . 

Finds out whether the intensity and area 
ratio of the match pair spot is very similar 
to each other. 
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SOLVE_ CLUSTER OPERATION 



Opens the match. cfg file to get the file name which will be 
processed. Opens all the .mmh, .cmp, .exm, .els, and .uni files. 

1. getpimw Fetches PI and MW values from input file. 

For each record in input file (.mmh), do the routines (2-5): 

2. get_record Gets the record from input file (.mmh). 

3. get__gid Finds out the gel index of a spot. 

If the spot is a unique spot, t hen call put record to write o ut 
the spot to output, otherwise, do the routines: check pattern 



4 . check_pattern 



Does a pattern matching using the current 
record. If it is not matched in the cluster 
pattern, call cluster_insert routine to 
insert this new cluster into the cluster 
list. If it is matched, but it is not exact 
match, call update_pattern routine to update 
the pattern of the cluster, and call 
record_insert routine to insert a record into 
a cluster. if it is matched and it is exact 
match, call record_insert routine to insert a 
record into a cluster. 



5. clusterjinsert Inserts a new cluster into the cluster list 



6. flush_pattern 



7 . sort cmp 



Outputs all the outstanding patterns/records. 
If the record in the cluster is an exact 
match, call exact_flush routine to output the 
record, otherwise call resolve_cluster 
routine to resolve the cluster and output the 
record. m case that it is a single 
reference match, don't output .els, .exm, snd 
. uni . 



Rearrange the order of 
increasing of PI value. 



, cmp file by 
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EXTRACT POTENTIAL MISMATCH 

Gets the matched records from .cmp. 

Reads one matched record from .smh each time. If 
there is no more record, the program is finished. 

Compares the record read in GetSmhRecord with 
record of .cmp by using the record number of the 
first gel (reference gel in single reference 
match) as key. If the comparison is same, call 
WritePcm, otherwise, call WriePmm. Jumps back to 
GetSmhRecord. 

Write the record read from .smh to .pern (Potential 
Correct Match) . 

Write the record read from .smh to .pmm (Potential 
Mis Match) . 
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* ad just, h : header file of ad just, c * 

* Author : Wenjeng Ko * 

* Date : November 1989 * 
************************^^ 

#define DOS 

#include <stdio.h> 
#include <ctype.h> 
#include <math.h> 

#ifdef DOS 

f include <stdlib.h> 

#include <alloc.h> 
#endif 

#def ine NEWLINE ' \ n ' 

#def ine EOS ' \0' 

#define BLANK ' r 



#define bool char 

#define TRUE l 

#def ine FALSE 0 

fdefine MAXLINE 128 

#define MAX_GELS 35 
fdefine MAXJ4ARKERS 50 

#define MK_BASE 9000 

#def ine X 0 

fdefine Y 1 

#define PI 0 

fdefine MW 1 

fdefine BELL 7 

fdefine SQUARE (x) ( ( x ) * (x) } 



struct { /* record #, pi, and mw of markers */ 

xnt rec; 

float pi,, mw; 
} mk_info [MAX_MARKERS] ; 



APPENDIX B-l 



WO 91/19274 



typedef struct tspot { 



PCT/US91/03620 



27 



xnt rec, 
Y/ 

newX, 
newY, 

mw, 

h, 

w; 

char name; 
float ii, 

area; 
struct tspot 
struct tspot 
SPOT; 



/* 
/* 
/* 
/* 
/* 



record number 
x-coord of the marker 
y-coord of the marker 
new x-coord of the spot 
new y-coord of the spot 
/* pi value 
/* mw value 
spot height 
spot width 
marker name 

intergrate intensity of the spot 
area 

/* pointer to next spot 
/* pointer to previous spot 



/* 

/* 

/* 

/* 

/* 
*next; 
*prev; 



*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 
*/ 



typedef struct tmarker { 
int rec, 

oldX, oldY, 
newX, newY, 
val, 
name ; 
float ii, 

area; 

struct tmarker *next; 
} L MARKER; 



/* record number of markers */ 



/* old X and Y values */ 

/* new X and Y values */ 

/* the value of the marker */ 

/* marker name */ 

/* intensity */ 

/* area ~ */ 

/* pointer to next marker */ 



typedef struct { 



rec, 


/* 


record number of the spot 


*/ 


oldX, 


/* 


old x-coord of the spot 


*/ 


oldY, 


/* 


old y-coord of the spot 


*/ 


newX, 


/* 


new x-coord of the spot 


*/ 


newY, 


/* 


new y-coord of the spot 


*/ 


ef f ect_r, 


/* 


effect radius 


*/ 


idx, idy, 


/* 


individual dx & dy of each marker 


*/ 


name; 


/* 


name 


*/ 



} MARKER; 



typedef struct gel { 
char name [10] ; 
int number_of_spots; 
SPOT *head; 

int x_low, x_high, y_low, y_high; 
int d_hpi , d_lpi , d_hmw , d_lmw# 
int d_hx , d_lx , d_hy , d_ly ; 
L__MARKER *Xmk_head; /* 
L_MARKER *Ymk_head; /* 
} GEL; 



/* gel name */ 
/* number of spots */ 
/* spot list header */ 
/* selected window */ 
/* delta of pi & mw */ 
/* delta of x & y */ 
X marker list header */ 
Y marker list header */ 
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* ad just. c : adjust the coordinate of each spot in study * 

* gel * * 

* Author : Wenjeng Ko * 

* Date : November 1989 * 
*********************** ***********************************/ 

#include "ad just. h" 

#include " d : \turboc\csslib\twindow . h" 



/*■ 
I 



subroutine protocol 



I 

-*/ 



void ReadGetCfgO ; 

void SetMarkerO; 

void Marklnsert () ; 

void AssignMarkers () ; 

void AdjustXYO; 

void PrintAttMks () / 

MARKER *GetNearestMark () ; 

MARKER *GetSecondNearestMark () ; 

L_MARKER *MarkAlloc ( ) ; 

SPOT *SpotAlloc () ; 



/*■ 
I 



global variables 



WINDOW 

FILE 

int 

GEL 

char 

int 

MARKER 



*msgwnd; /* message window pointer 

*Fp_com; /* file pointer for common use 

Num_gels =0; /* number of gels 

Gels[MAX_GELS] ;/* gel information 

File_dir[40] ; /* file directory 

Num_mk; /* number of markers 

Markers [MAX_GELS] [MAX_MARKERS] ; /* marker info, 



.*/ 

*/ 

*/ 

*/ 

*/ 

*/ 
*/ 

*/ 
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/* */ 

main (argc, argv) 
int argc; 
char *argv[]; 
{ 

int i, j; 

char cmd[80], cfg[55]; 

char gelname[13], lstfile[55], mksfile[55], attfile[55]; 
char file_dir [40] , str[80]; 

if (argc != 2) error ("Usage : ", "adjust data_directory") ; 
strcpy (File_dir , argv [ 1 ] ) ; 
Num_gels = ReadMatchCf g () ; 
ReadGetCfg () ; 

for (i « 0; i < Num_gels; i++) 

Gels [i] .number_of_spots = ReadListFiles (i) / 

/* Adjust x and y coordinates according to Ref. gel (gel[0]) */ 
AssignMarkers (0) ; 
SetMarkersMove (0) ; 
SetRefNewCoor () ; 

for (i = 1; i < Num_gels; ( 
AssignMarkers (i) ; 
SetMarkersMove (i) ; 
AdjustXY(i) ; 

> 

for <i = 0; i < Num_gels; i++) { 
AssignMkValues (i) ; 
AssignNewXy (i) ; 

} 

/* Get the .mks and .att files in this step */ 
for (i = 0; i < Num_gels; { 
GetPiMw(i) ; 

} 

/* Print, out the .att and .mks files */ 
for (i = 0; i < Num_gels; i++) { 
Print At tMks (i) ; 

} 

} /* main */ 
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/* 

AssignNewXy (idx) 
int idx; 
{ 



30 



*/ 



int 

L MARKER 



if j; 

*mp; 



mp = Gels [idx] .Xmk_head; 
while (mp != NULL) { 

for (i = 0; i < Num_mk; i++) { 

if (mp->rec ~ Markers [idx] [i] *rec) { 
mp->newX = Markers [idx] [i] .newX; 
mp->newY = Markers [idx] [i] .newY; 
break; 

} 

} 

mp = mp->next; 
} /* while */ 



mp « Gels [idx] .Xmk_head; 
j - 0; 

while (mp != NULL) { 
if (j — 0) { 

Gels [idx] .d_hx = mp->newx 
Gels [idx] .djtipi = mp->val 

. > . 

else if (j — Num_mk -2) { 
Gels [idx] .d_lx = mp->newX 
Gels [idx] .d_lpi « mp->val 

} 



mp->next->newX; 
mp->next->val; 



mp->next->newX; 
mp->next ->val ; 



} 



mp = mp->next; 
/* while */ 



mp = Gels [idx] .Ymk_head; 
while (mp != NULL) { 

for (i = 0; i < Num_mk; i++) { 

if (mp->rec — Markers [icix] [i] .rec) { 
mp->newX = Markers [idx] [i] .newX; 
mp->newY - Markers [idx] [i] .newY; 
break; 

1 



} 



} 

mp = mp->next; 
/* while */ 



mp - Gels [idx] .Ymk_head; 
j = 0; 

while (mp != NULL) { 
if (j == 0) { 
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Gels [idx] .d_hy « mp->newY - mp^>next->newY; 
Gels [idx] .d_hmw = mp->next->val - mp->val; 

} 

else if (j — Num_mk -2) { 

Gels [idx] .d_ly = mp->newY - mp->next->newY; 
Gels [idx] .dJLmw = mp->next->val - mp->val; 

} 

mp = mp->next; 

} /* while */ 
} /* AssignNewXy */ 



/* . . */ 

/* 

ReadMatchCfg Open match. cfg and get the Gel__list 

(gel name table) . 

*/ 

static int ReadMatchCfg () 
{ 

char line [MAXLINE+1] , dm [15]; 
int i, num; 



/* Get the gel name from match, cfg file */ 
tifdef UNIX 

sprintf (line, "%s/ /match. cfg", File_dir) ; 
#endif 
#ifdef DOS 

sprintf (line, ,, %s\\match.cfg", File_dir) ; 
#endif 

if ((Fp_cora = fopen(line, "r") ) == NULL) 
error ("Can' t open file match. cfg "!"); 

for ( i - 0/ i < 4/ i++) 

fgets (line, MAXLINE, Fp_com) / 
num =0; 

while (f scanf (Fp_com, "%s n , Gels [num] .name) != EOF) { 
if (strlen (Gels [num] .name) > 6) 

error ("Gel name is over 6 characters — > 
Gels [num] .name) ; 

num++; 

} 

f close (Fp_com) / 
return num; 
} /* ReadMatchCfg */ 



/* */ 

/* 

ReadGetCfg read marker infomation ( get. cfg ) 
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*/ 

void ReadGetCfgO 
{ 

L_MARKER *mtp; 
char str [MAXLINE] ; 
float pi, mw; 

int i, j, chkcnt, rec, nummk; 

#ifdef UNIX 

sprintf (str, "%s//get . cfg" , File dir) / 
#endif 
#ifdef DOS 

sprintf (str, "%s\\get . cfg", File dir> ; 
#endif "~ 

if ( (f reopen (str, fl r", Fp_com) ) — NULL) 
error ("Can' t open file str); 

f gets (str, MAXLINE, Fp_com) ; 
sscanf(str, "%d ", & nummk ) ; 

for (i = 0; i < nummk; i++) { 
f gets (str, MAXLINE, Fp_ com) ; 
sscanf(str, »%d %f %f », &rec, &pi, &mw) ; 
mk_info[i] .rec = rec; 
mk_info[i] .pi « pi; 
mk_inf o [ i ] . mw = mw; 

} 

f close (Fp_com) ; 
} /* ReadGetCfg */ 



ReadListFiles read spot information from .1st file 

*/ 

int ReadListFiles (idx) 

int idx; 

{ 

SPOT *stemp; 

char imagename[15] , str [MAXLINE] ; 
int i, rec, x, y, spotname, h, w, num; 
float ii, area; 
bool is_marker; 

#ifdef UNIX 

sprintf (str, "%s//%s . 1st" , File_dir, Gels [idx] .name) ; 
#endif 
#ifdef DOS 
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sprintf (str, "%s\\%s . 1st", File_dir, Gels[idx] .name) ; 
#endif 

if ( (f reopen (str, "r", Fp_com) ) == NULL) 
error ("Can't open file str); 

ford = 1; i <= 5; i++) fgets(str, MAXLINE-1, Fp_com) ; 
• num = 0 / 

while ((f gets (str, MAXLINE, Fp_com) ) 0) { 
is_marker = FALSE; 
if (Str [5] == BLANK) 

sscanf(str, f, %d %s %d %d %f %f %d %d", Srec, imagename, 

&x, &y, &ii, &area f &h, &w) ; 

else { 

is_marker = TRUE/ 

sscanf (str, "%d %d %s %d %d %f %f %d %d", &rec, Sspotname, 
imagename, &x, &y, &ii, &area, &h, &w) ; 

} /* else */ 
num++; 

if (is_marker — TRUE) rec = spothame + MK__BASE; 

stemp - SpotAllocO; 

stemp->rec = rec; 

stemp->x = x; 

stemp->y = y; 

stemp->ii - ii; 

stemp->area = area; 

stemp->h = h; 

stemp->w = w; 

Spotlnsert (idx, stemp) ; 

if <is_marker == TRUE) SetMarker (idx, stemp); 
} /* while */ 
f close (Fp_com) ; 
return num; 
} /* ReadListFiles */ 



/*.... */ 

/* 

Spotlnsert insert the node to list in decreasing order 

This is a double link list. 

*/ 

static Spotlnsert (which_gel, stemp) 
int whi ch_ge 1 ; 
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SPOT *sterop; 
{ 

SPOT *tp, *head; 

head = Gels [which_gel] .head; 

if (head == NULL) { 
stemp->prev «= NULL; 
Gels [which_gel] .head = stemp; 
return; 

> 

tp = head; 

while (tp ! = NULL) { 

if (stemp->x <= tp->x) { 

if (tp->next — NULL) { /* this is the last item */ 
stemp->prev = tp; 
tp->next = stemp; 
break; 

1 

else { /* this is not the last item */ 

if (stemp->x >= tp->next->x) { 
stemp->prev = tp; 
stemp->next = tp->next; 
tp->next->prev =* stemp; 
tp->next = stemp; 
break; 

} 

else { /* go to next loop */ 
tp = tp->next; 

} 

} 

} /* if */ 

else { /* put in the first position */ 
stemp->next = tp; 
stemp->prev = NULL; 
tp->prev = stemp; 
head = stemp; 
break; 

} 

} /* while */ 

Gels [which_gel] .head = head; 
} . /* Spot Insert */ 

/* */ 

/* 

SetMarker link the marker to marker list 

*/ 
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void SetMarker (which_gel, stemp) 
int which_gel; 
SPOT * stemp ; 
{ 

L_MARKER *mtemp; 

mtemp = MarkAllocO; 
mtemp->rec = stemp->rec; 
mtemp->oldX = stemp->x; 
mtemp->oldY = stemp->y; 
mtemp->val = stemp->pi; 
mtemp->name = stemp->name; 
Marklnsert (which — gel, mtemp, PI) ; 

mtemp = MarkAllocO; 
mtemp->rec = stemp->rec; 
mtemp->oldX = stemp->x; 
mtemp->oldY = stemp->y; 
mtemp->val = stemp->mw; 
mtemp->name = stemp->name; 
Marklnsert (which_gel, mtemp, MW) ; 
} /* SetMarker */ 



/* 
/* 

Marklnsert insert marker to marker list 

*/ 

void Marklnsert (whichjgel, mtemp, type) 

LJ4ARKER *mtemp; 

int type, which_gel; 

i 

L_MARKER *head, *tp, *prev; 

if (type == PI) head = Gels [which_gel] .Xmk_head; 

else head = Gels [which_gel] . Ymk_head; 

tp = head; 

if (tp == NULL) head = mtemp; 
else { 

prev = head; 

while (tp != NULL) { 

if (type == PI tp->oldX < mtemp->oldX) break; 
if (type == MW && tp->oldY < mtemp->oldY) break; 
prev = tp; 
tp = tp->next; 
} /* while */ 

if (tp == prev) { 
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mtemp->next = tp; 

head = mtemp; 
} /* if */ 
else { 

mtemp->next = tp; 

prev->next = mtemp; 
} /* else */ 
} /* else */ 

if (type « PI) Gels [which_gel] .Xmk_head - head; 
else Gels [which_gel] .Ymk_head ■ head; 

} /* Marklnsert */ 

/* */ 

void AssignMarkers (which_gel) 

int which_gel; 

{ 

L_MARKER *m; 

int xlow, xhigh, ylow, yhigh, mk_number; 

static int first = TRUE; 

xlow «= ylow = 9999; 

xhigh = yhigh = -1; 

m = Gels [which_gel] .Xmk_head; 

mk__number = 0; 

while (m NULL) { 

if (m->oldX < xlow) xlow = m->oldX; 

if (m->oldX > xhigh) xhigh = m->oldX; 

if (m->oldY < ylow) ylow = m->oldY; 

if (m->oldY > yhigh) yhigh = ra->oldY; 

Markers [which_gel] [mk_number] .rec = m->rec; 
■ Markers [which_gel] [mk_number ] .oldX = m->oldX; 
Markers [which_gel] [mkjiumber] ,oldY = m->oldY; 
Markers [which_gel] [mkjiumber] ,1131116 = Nuinjtik; 

mk_number ++ ; 
m = m->next ; 

} 

if (mk_number < 3) 

error ("Number of markers should be over 3« # " ."); 
if (mk_number > MAX__MARKERS ) 

error ( "Number of markers should be less than 50", "-"); 
if (first) { 

Num_mk = mk_number; 

first = FALSE; 

} 

else { 

if (rak_number != Num_mk) 

error (Gels [which_gel] .name, 
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" has different marker number."); 

) 

Gels [which_gel] .x_low = xlow; 
Gels [which_gel] .x_high = xhigh; 
Gels [which_gel] .y_low = ylow; 
Gels [which_gel] .y_high = yhigh; 
'/* AssignMarkers */ 
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/* • */ 

/* ' 

SetMarkersMove : Compare markers in Gel i with Gel 0, 

put the difference in idx and idy. 

SetMarkersMove (which_gel) 

int which_gel; 

{ 

MARKER m; 
int i, j; 

if (which_gel == 0) { /* This is ref gel */ 

for (i = 0; i < Num_mk; i++) { 

Markers [0] [i] .newX = Markers [0] [i] .oldX; 
Markers [0] [i] .newY - Markers [0] [i] .oldY; 
Markers [0] [i] .idx = Markers [0] [i] .idy - 0; 

return; 

> 

for (i - 0; i < Num_mk; i++) { 
• for (j = 0; j < Num_mk; j++) { 

if (Markers [which_gel] [i] .rec == Markers[0] M] .reel 
break; 

} 

Markers [which_gel] [i] .idx = 

Markers [0] [j] .oldX - Markers [which gel] [i] .oldX; 
Markers [which_gel] [i] .idy = 

Markers [0] [j] .oldY - Markers [which_gel] [i] .oldY; 
Markers [which_gel] [i] .newX = Markers [0] [ j] .oldX; 
^ Markers [which_gel] [i] .newY = Markers [0] [j] .oldY; 

} /* SetMarkersMove */ 



void AdjustXY(which_gel) 

int which gel; 

{ 

SPOT *sp; 

sp = Gels [which_gel] .head; 
while (sp != NULL) { 

if (sp->x <= Gels[which_gel] .x high && 
sp->x >= Gels[which_gel] .x__low && 
sp->y <= Gels[which_gel] .y_high &'& 
sp->y >= Gels[which_gel] .y__low ) 
Adjust (which_gel / sp) ; 
sp = sp->next; 

} 
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} /* AdjustXY */ 

/* */ 

/* 

SetRefNewCoor Set the reference gel's new X, Y value. 

*/ 

SetRefNewCoor () 
{ 

SPOT *sp; 

sp - Gels [0] .head; 
while (sp != NULL) { 

if <sp->x <= Gels [0] .x_high && 
sp->x >= Gels [0] .x_low && 
sp->y <= Gels [0] .y_ high && 
sp->y >= Gels [0] .y_JLow ) { 
sp->newX = sp->x; 
sp->newY = sp->y; 

> 

sp = sp->next; 

} 

} /* SetRefNewCoor */ 



/* */ 

/* 

Adjust Adjust the spots in the study gels to it's markers 

*/ 

static Adjust (which_gel, sp) 
int whi ch_ge 1 / 

SPOT *sp; 
{ 

MARKER *mark0, *markl, *mark2; 

int tl, t2, t3, t4, tx, ty; 

markO = GetNearestMark (which_gel, sp) ; 
markl = GetSecondNearestMark (which_gel, sp) ; 

tl = mark0->oldX - sp->x; 
t2 = mark0->oldY - sp->y; . 
t3 — markl->oldX - sp->x; 
t4 = markl->oldY - sp->y; 

if ( SQUARE (tl) +SQUARE (t2) <= SQUARE ( (long) mark0->ef f ect_r) ) { 
tx = mark0->idx; 
ty = mark0->idy; 

} 

else if ( (SQUARE (tl) + SQUARE (t2) <= 4 * 
SQUARE ( (long) mark0->ef f ect_r ) ) 
(SQUARE (t3) + SQUARE (t4) <= 4 * 
SQUARE ( (long) markl->ef f ect_r) ) ) { 
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+ markl->idx) / 2.0; 
+ markl->idy) / 2.0; 



else { 

tx = mark0->idx; 
ty = mark0->idy; 

} 



sp->newX = sp->x + 
sp->newY = sp->y + 



tx; 
ty; 



} /* Adjust */ 



/* 



*/ 



MARKER 

int 

SPOT 



*GetNearestMark(which_gel, node) 

which_gel; 

*node; 



int i; 

unsigned long min = 9999999; 
long dist; 
long tpl, tp2; 
int mark; 

for ( i = 0 ; i < Numjnk ; i++ ) { 

/* Since all Study-gels' marker. newX and marker. newY are 
equal to Ref-gel's, and we use newX, newY to find the 
nearest marker, therefore gel [Ref_gel] ' s marker list 
can be use for every Study-gel's marker list 

*/ 

tpl * node->x - Markers [which_gel] [i] .oldX; 
tp2 = node->y - Markers [which_gel] [i] .oldY; 

dist - SQUARE (tpl) + SQUARE <tp2>; 

if ( dist < min ) { 
min = dist; 
mark - i; 



} 



return ( ^Markers [which_gel] [mark] )j 
} /* GetNearestMark */ 



/* 



*/ 



MARKER *GetSecondNearestMark(which_gel, node) 

int which_gel; 

SPOT *node; 
{ 
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int i; 

unsigned long mini = 9999999, min2 = 9999999/ 
long dist; 
long tpl, tp2; 
int markl, mark2; 



for ( i = 0 ; i < Num_mk ; i++ ) { 

/* Since all Study-gels' marker. newX and marker. newY are 
equal to Ref-gel's, and we use newX, newY to find the 
nearest marker, therefore gel [Ref_gel] ' s marker list 
can be use for every Study-gel's marker list 

*/ 

tpl = node->x - Markers [which_gel ] [i] . oldX; 
tp2 = node->y - Markers [which_ge ; ] [i] .oldY; 

dist « SQUARE (tpl) + SQUARE (tp2); 

if ( dist < mini ) { 
min2 = mini; 
mark2 = markl; 
mini = dist; 
markl = i; 

} 

else if (dist < min2) { 
min2 = dist; 
mark2 = i; 

} 

} 

return ( ^Markers [which_gel] [mark2] ) ; 
} /* GetSecondNearestMark */ 



GetPiMw get PI & MW values for each spot 

*/ 

GetPiMw (which_gel) 
int which gel; 

{ 

SPOT *sp; 



sp = Gels [which__gel] .head; 
while ( sp. NULL ) { 

if( sp->x >= Gels[which_gelJ f x low && 
sp->x <= Gels [which_gel} .x~high 
sp->y >= Gels [which_gel] .y_low && 
sp->y <= Gels [which_gel] .y_high ) { 
CalPi (which_gel, sp) ; 
CalMw ( whi ch_ ge 1 , sp ) ; 
} /* if */ 
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else { 

sp->pi = -1000; 
sp->mw = -1000/ 

} 

sp = sp->next; 
} /* while */ 
} /* GetPiMw */ 



/* 

/* 

CalPi calculate the pi value of the spot 

*/ 

CalPi (which_gel, sp) 
int whi ch_ge 1 / 

SPOT *sp; 
{ 

L__MARKER *mp; 
int i; 
long t; 



mp = Gels [which gel] .Xmk head; 
i = 0/ 

while <mp != NULL) { 

if (sp->newX>mp->newX && i == 0) { /* over right margin */ 
t = ( (long) sp->newX - mp->newX) * 
Gels[which_gel] .d_hpi / 
Gels [which_gel] .d_hx - mp->val; 
sp->pi = (int) t; 
break; 

} 

else if (sp->newX == mp->newX) { 
sp->pi = mp->val; 
break; 

} 

else { 

if (mp->next != NULL) { 

if (sp->newX > mp->next->newX) { 

t = ( (long) sp->newX - mp->next->newX) * 
(mp->val - mp->next->val) / 
(mp->newX - mp->next->newX) ; 
sp->pi « mp->next->val + (int) t; 
break; 
} ■ 

} 

else { 

t = ((long) k. ->newX - sp->newX) * 

Gels [whio:s_gel] .d_lpi / 

Gels[whi. :_gel] .d_lx; 
sp->pi = mp->v=il - (int) t; 
break; 
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} 

} 

mp = mp->next; 

i++; 
} /* while */ 
/* CalPi */ 
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'* */ 

/* 

CalMw calculate the mw value of the spot 

*/ 

CalMw (which_gel, sp) 
int which__gel ; 

SPOT *sp; 
{ 

L_MARKER *mp; 
int i; 
long t; 



mp = Gels [which gel] .Ymk head; 
i = 0; - 
while (mp != NULL) { 
if (sp->newY>mp->newY i — 0) { /* over right margin */ 
t = ( (long) sp->newY - mp->newY) * 
Gels [which_gel] .d_hmw / 
Gels[which_gel] .d_hy; 
sp->mw = mp->val - (int) t; 
break; 

} 

else if (sp->newY — mp->newY) { 
sp->mw = mp->val; 
break; 

} 

else { 

if (mp->next != NULL) { 

if (sp->newY > mp->next->newY) ( 

t = ( (long) sp->newY - mp->next->newY) * 
(mp->val - mp->next->val) / 
(mp->newY - mp->next->newY) ; 
sp->mw « mp->next->val + (int) t; 
break; 

} 

} 

else { 

t = ( (long) mp->newY - sp->newY) * 

Gels[which_gel] .d_lmw / 

Gels [which_gel] .d_ly; 
sp->mw = mp->val - (int) t; 
break; 

} 

} 

mp = mp->next; 
i++; 
} /* while */ 
} /* CalMw */ 
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/* */ 

/* 

AssignMkValue — assign the value of marker to PI marker list 

*/ 

static AssignMkValues (which_gel) 

int which_gel; 

{ 

L_MARKER *mp; 

/* Assign pi values */ 

mp = Gels [which_gel] .Xmk_head; 

while (mp != NULL) { 

mp->val = (int) (mk_inf o [mp->rec-MK_BASE-l] .pi * 100); 

mp = mp->next; 

} 

/* Assign mw values */ 

mp = Gels [which_gel] . Ymk_head; 

while (mp != NULL) { 

mp->val = (int) (mk_info [mp->ree-MK_BASE-l] .mw * 100); 

mp = mp->next; 

} 

} /* AssignMkValues */ 



/* */ 

/* 

MarkAlloc allocate the marker node in memory 

*/ 

Static L_MARKER *MarkAlloc() 
{ 

L_MARKER *mtemp; 

if< (mtemp = (L_MARKER *) malloc (sizeof (L_MARKER) ) ) == NULL ) 
error ("Out of memory","!"); 

mtemp->val = -1000; 
mtemp->next = NULL; 
return (mtemp) ; 
) /* MarkAlloc */ 



/* . */ 

/* 

SpotAlloc allocate the spot node in memory 

*/ " 

SPOT *SpotAlloc() 

{ 

SPOT *stemp; 
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if( (stemp = (SPOT *) malloc (sizeof (SPOT) ) ) 

error ( "Out of memory" ,"!"); 
stemp->pi = -1000; 
stemp->mw = -1000/ 
stemp->newX = -1; 
stemp->newY = -1; 
stemp->next = NULL; 
stemp- >prev = NULL; 

return (stemp) ; 
} /* Spot Alloc */ 



/* */ 

void Print At tMks (which_gel) 

int. which_gel; 

{ 

SPOT * stemp , *stp; 

LJMARKER *mtemp, *mtp, *tmp [MAX_MARKERS] ; 

int i; 

char str[55]; 

/* Print out the .att files */ 
#ifdef UNIX 

sprintf(str, "%s//%s . att" , File_dir, Gels [whichjgel] .name) ; 
fendif 
#ifdef DOS 

sprintf(str, "%s\\%s.att", File_dir, Gels [which gel] .name) ; 
#endif ~ 

if ( (f reopen (str, "w", Fp_com) ) — NULL) 
error ("Can't open file ", str) ; 

fprintf (Fp_com, "Filename: %s .att\n\n", Gels [which_gel] .name) ; 
fprintf (Fp_com, 

"Rec# X Y NewX NewY PI M W II Area Ht 

Wd_\n") ; ~ — 

stemp = Gels [which_gel] .head; 
while ( stemp NULL ) { 

if ( stemp->pi != -1000 stemp->mw != -1Q00 ) { 
fprintf (Fp_com, 

"%4d %4d %4d %4d %4d %6.2f %6.2f %7.3,f %7.3f %4d %4d\n" / 
stemp->rec, stemp->x # stemp->y # 
stemp->newX, stemp->newY, 

stemp->pi/100 .0, stemp->mw/100 .0, stemp->ii, 
stemp->area, stemp->h, stemp->w) ; 

> 

stp = stemp; 

stemp = stemp->next; 
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} /* while */ 
f close (Fp_com) ; 

/* Print out the .als (modified .1st) files. The spots outside 

window will be ignored. 
*/ 

#ifdef UNIX 

sprintf(str, "%s//%s . als" , File_dir, Gels [which_gel] .name) / 
#endif 
#ifdef DOS 

sprintf (str, "%s\\%s .als", File_dir, Gels [which_gel] .name); 
#endif 

if ( (f reopen (str, "w", Fp_com) ) == NULL) 
error ("Can't open file ", str) / 

fprintf (Fp_com, "Filename: %s.als\n", Gels [which_gel] .name) ; 
fprintf (Fp_com, "Window = X : %4d - %4d\n" / 

Gels [which_gel] .x_low, Gels [which_gel] .x_high) ; 
fprintf (Fp_com, " Y : %4d - %4d\n\n", 

Gels [which_gelj .y_low, Gels [which_gel] .y_high) / 
fprintf (Fp_com, 

"Rec# Spotname Image NewX NewY II Area Ht 

Wd \n") ; 

stemp = Gels [which_gel] .head; 
while ( stemp != NULL ) { 

if( stemp->pi ! = -1000 && stemp->mw ! = -1000 ) { 
if (stemp->rec >= MKJBASE) 
. . fprintf (Fp_com, "%4d %-10d %-10s %4d %4d %7.3f %7.3f %4d 

%4d\n", 

stemp->rec, stemp->rec-MK_BASE, 
Gels [which_gel] .name, 
stemp->newX, stemp->newY, stemp->ii, 
stemp->area, stemp->h, stemp->w) ; 

else 

fprintf (Fp_com, "%4d %-10s %-10s %4d %4d %7.3f %7.3f %4d 

%4d\n", 

stemp->rec, " V, Gels [which_gel] .name, 
stemp->newX, stemp->newY, stemp-r^ii, ' 
stemp->area, stemp->h, stemp->w) ; 

} • 

stp = stemp; 

stemp = stemp->next; 

free (stp) ; 
} /* while */ 
f close (Fp_com) ; 

/* Print out the .mks files */ 
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#ifdef UNIX 

sprintf (str, "%s//%s .inks" , File_dir, Gels [which gel] .name) ; 
#endif "* 
frifdef DOS 

sprintf (str, n %s\\%s .niks", File dir, Gels [which gel] .name) ; 
#endif " - 

if ( (f reopen (str, "w", Fp_com) ) == NULL) 
error ("Can 7 t open file str); 

fprintf (Fp_com, "%s .mks\n n , Gels [which_gel] .name) ; 
f printf (Fp_com, "Window = X : %4d ~ %4d\n", 

Gels [which_gel] .x_low, Gels [which_gel] .x_high) ; 
fprintf (Fp_com, " Y : %4d - %4d\n ,f , 

Gels [which_gel] .y_low, Gels [which_gel] .y_high) ; 

/* sort the markers by record number */ 
mtemp - Gels [which_gel] .Xmk_head; 
for( i = 0; i < Num_mk; i++ ) { 

tmp [mtemp->rec - MK_BASE] = mtemp; 

mtemp = mtemp->next; 
} /* for */ 

fprintf (Fp_com, 

"\n %2d Markers \n", Numjnk) ; 

fprintf (Fp_com, "Marker_name X Y \n") ; 

for( i ■ 1; i <= Num_mk; i++ ) { 

fprintf (Fp_com, « %4d %4d %4d\n M , 

tmp[i]->rec, tmp [i] ->oldX, tmp [i] ->oldY) ; 
} /* for */ 
f close (Fp_com) / 
} /* PrintAttMks */ 



error give the error message 

*/ 

error (si, s2) 
char *sl, *s2; 
{ 

error_message(sl, s2, 0, OK); 
exit (1) ; . 
} /* error */ 

/*==""— ==——==== s end of ad just, c ==================* / 
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Documentation Conventions 

.1st = raw data file (X-Y coordinate) 

.als = adjusted data file (X-Y coordinate) after run "match" 

.att = adjusted data file (PI-MW coordinate) after run "match" 

.mks = information about the markers 

•mmh = multiple reference matched image file 

.smh = single reference matched image file 

.cmp = composite match image file 

.pno « reference file for reading .cmp map 

.exm = exact match spots file 

.els = cluster spots file 

.uni = unique spots file 

.pmm » potential mismatched spots 

.pern = potential correct matched spots 

.lmt « metafile for drawing spot list map (X-Y coordinate) 
. .amt = metafile for drawing spot list map (PI-MW coordinate) 
•cmt = metafile for drawing compositive map 

.mmt = metafile for drawing the moving vector in any two 
matched gels 

xdrawout.plt « output file for HPGL and PostScript of printer 
or plotter 

get.cfg = configuration file of get pots 
match. cfg = configuration file of match 

draw.cfg - configuration file of comparison of .att and .cmp 
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559 
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0.601 
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670 
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841 
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204 
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244 
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93 

75 
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52 
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53 
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54 

94 
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77 

184 
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Print file name: a.att, Time*: 15:12 n a te: 17^/1 oao 



Filename : a . att. 



PCT/US91/03620 



-Page ; 



Rec# 

9001 
136 
150 
113 
137 
112 
110 
157 
149 
73 
107 
135 
109 
148 
108 
105 
134 

9002 
106 

9003 
133 
146 
155 
104 
132 
102 
145 
103 
131 
100 

9004 
101 
97 
128 
144 
98 
99 
129 
127 
96 
95 

9005 

9006 
94 



696 
696 
694 
690 
658 
650 
636 
636 
630 
628 
620 
620 
600 
586 
576 
570 
560 
556 
524 
512 
504 
500 
492 
490 
480 
478 
460 
454 
452 
432 
426 
416 
380 
376 
372 
366 
334 
332 
324 
322 
292 
264 
240 
240 



484 
357 
430 
319 
369 
27 4 
308 
457 
402 
250 
268 
348 
316 
389 
302 
265 
336 
447 
302 
252 
372 
417 
464 
319 
333 
277 
417 
298 
324 
261 
360 
317 
259 
349 
385 
294 
291 
382 
332 
261 
290 
250 
346 
312 



NewX 
696 
696 
694 
690 
658 
650 
636 
636 
630 
628 
620 
620 
600 
586 
576 
570 
560 
556 
524 
512 
504 
500 
492 
490 
480 
478 
460 
454 
452 
432 
42 6 
416 
380 
376 
372 
366 
334 
332 
324 
322 
292 
264 
240 
240 



NewY 
484 
357 
430 
319 
369 
274 
308 
457 
402 
250 
268 
348 
316 
389 
302 
265 
336 
447 
302 
252 
372 
417 
464 
319 
333 
277 
417 
298 
324 
261 
360 
317 
259 
349 
385 
294 
291 
382 
332 
261 
290 
250 
346 
312 



-1.00 
-1.00 
-1.02 
-1.05 
-1.28 
-1.33 
-1.43 
-1.43 
-1.48 
-1.49 
-1.55 
-1.55 
-1.69 
-1.79 
-1.86 
-1.90 
-1.98 
-2.00 
-2.18 
-2.24 
-2.32 
-2.35 
-2.42 
-2.44 
-2.53 
-2.55 
-2.70 
-2.76 
-2.78 
-2.95 
-3.00 
-3.07 
-3.29 
-3.31 
-3.34 
-3.38 
-3.57 
-3.59 
-3.63 
-3.65 
-3.83 
-4.00 
-4.21 
-4.21 



19.00 
21.76 
20.61 
24.20 
21.39 
26.68 
24.81 
20.02 
20.97 
28.00 
27.01 
22.54 
24.37 
21.13 
25.14 
27.18 
23.27 
20.39 
25.14 
27 . 8 9 
21.35 
20.78 
19.76 
24.20 
23.43 
26.52 
20.78 
25.36 
23.93 
27.40 
21.50 
24.31 
27.51 
22.46 
21.19 
25.58 
25.75 
21.22 
23.49 
27.40 
25.80 
28.00 
22.71 
24.59 



II 



3. 
0. 
1. 
3. 
0. 
1. 
2. 
1. 
1. 
1. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
4. 
1. 
8. 
3. 
0. 
0. 
1. 
0. 
1. 
0. 
0. 

1. 

0. 
4. 
0. 
0. 
3 
0 
0 
0 
0 
0 
0 
0 
3 
7 
0 



652 
300 
023 
868 
859 
391 
391 
162 
.845 
.006 
.049 
.281 
.022 
.264 
.111 
.052 
.055 
.757 
.773 
.928 
.676 
.203 
.379 
.404 
.039 
.250 
.140 
.437 
.888 
.884 
.162 
.131 
.702 
.781 
.051 
.026 
.162 
.074 
.191 
.150 
.296 
.127 
.771 
.531 



Area 

5.140 



Ht Wd 



1. 
2. 
5. 
1. 
1, 



1. 
1, 
0, 

0, 
5, 



.030 
.620 
.060 
.990 
.670 
3.340 
2.200 
2.860 
1.800 
0.190 
1.140 
0.190 
.540 
.090 
.720 
.480 
.120 
2.330 
6.710 
5.800 
1.300 
1.910 
1.380 
0.080 
1.540 
1.140 
0.980 
2.170 
1.910 
3.900 
0.610 
1.960 
5.190 
0.370 
0.270 
1.110 
0.930 
1.300 
0.980 
1.540 
4.590 
7.790 
1.190 



36 
10 
22 
32 
20 
22 
24 
28 
24 
20 

4 
14 

4 
18 
12 
10 

8 
32 
20 
38 
36 
14 
16 
14 

2 
20 
14 
12 
20 
18 
32 
10 
20 
30 

6 

4 
14 
12 
12 
16 
16 
26 
32 

8 



28 
18 
30 
36 
20 
22 
30 
24 
26 
20 

8 
16 

8 
16 
18 
14 
12 
32 
20 
46 
38 
18 
32 
20 

6 
16 
16 
16 
20 
18 
28 
12 
20 
32 
12 
10 
14 
14 
18 
14 
18 
38 
52 
26 
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Print file name: a.als, Time: 15:19 DatP: 12/1/1989 gage: 

Filename: a.als 
Window =-X : 240 - 696 
Y : 250 ~ 484 



Rec# 


Spot name 


Image 


NewX 


NewY 


II 


Area 


Ht 


Wrf 

"U. 


9001 


1 
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484 


3 . 652 


5 .140 


~36 


28 
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696 


357 


0 .300 
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10 


18 


150 
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694 
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30 
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690 
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•J X 


3 fi 
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0 . 859 
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9 0 
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9 0, 
x U 
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1 670 
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JO 


157 




d 


636 


4 57 


1 1 69 
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9 9 A A 
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24 
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v J U 


*» V £. 


1 Q45 
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73 






628 


2 50 




1 Q A A 
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-J X V 
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0 9 64 


x a 34U 


"t D 
1 O 


16 


108 
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57 6 
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0 111 

v • X X X 


1 AQft 

x . u y u 


1 9 


1 O 
1 0 


105 




a 


570 


265 


0 0 52 


0 79 n 

u . / z u 


xU 


14 


134 




a 




336 


v • UD3 


A /ion 
V . 4oU 


Q 
O 


12 


9002 


2 


a 


33 O 


44*7 


•1 • / D / 


c n oa 
3 • 1^0 


32 


32 


106 






59 £ 


3 f!9 
O v ^ 


"1 9*7 ^ 
x . / / J 




20 


20 


9003 


3 


a 


512 


252 


8 92ft 


0 . / x U 


^ ft 
JO 


A C 
4 O 


133 




a 


504 


372 


3 676 


5 800 


3 6 
O O 


"3 ft 
J O 


146 




a 


500 


417 


0 203 


1 300 

1 .OUU 


1 4 
x 4 


1 ft 


155 




a 


492 


464 

1 V? *1 


0 37 9 


1 Q1 A 


Id 


9 9 
JZ 


104 




B 




31 Q 


1 AHA 
x . 4 U 4 


t *a q n 
1 . Jo U 


14 


20 


132 




a 
<a 




333 




A Ck O A 


2 


6 


102 




3 


4 / O 


9*7*7 


T 9 c A 


1 . 540 


20 


16 


145 

-i- I -J 




a 


46 n 

40U 


4X / 


A 1 /I A 


1 . 140 


14 


16 


-i. V J 




a 


A ^ A 
404 


9 a q 


0 - 437 


0 . 980 


12 


16 


X .J X 




2 




J^4 


1.888 


2 . 170 


20 


20 


100 




d 


4 3 9 


401 


A O D /I 

U . bo 4 


1 . 910 


18 


18 


9004 


4 


a 


426 


3 60 


A 1 £9 
4 . x Ox 


*5 a a a 

j . you 


J2 


28 


101 




a 


416 


317 


0 . 131 


0 610 
v • o x \j 


1 0 
X u 


1 9 
xx 


97 




a 


380 


259 


0.702 


1.960 


20 


20 


128 




a 


376 


349 


3.781 


5.190 


30 


32 


144 




a 


372 


385 


0.051 


0.370 


6 


12 


98 




a 


366 


294 


0.026 


0.270 


4 


10 


99 




a 


334 


291 


0.162 


1.110 


14 


14 


129 




a 


332 


382 


0.074 


0.930 


12 


14 


127 




a 


324 


332 


0.191 


1.300 


12 


18 


96 




a 


322 


261 


0.150 


0,980 


16 


14 


95 




a 


292 


290 


0.296 


1.540 


16 


18 


9005 


5 


a 


264 


250 


3.127 


4.590 


26 


38 


9006 


6 


a 


240 


346 


7.771 


7,790 


32 


52 


94 




a 


240 


312 


0.531 


1.190 


8 


26 
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Py;lnt file name: h.atit. £jjnej 15-14 n^ e! i^V/nopo 
Filename: b.att 



Rec# 

9001 
178 
198 
175 
213 
220 
176 
212 
211 
139 
195 
140 
174 
197 
173 
196 
137 

9002 
170 
171 
172 

9003 
193 
208 
194 
209 
168 
191 
218 
166 
207 
135 
169 
192 
167 
189 
134 

9004 
165 
188 
186 
131 
163 
206 
130 
162 
187 
164 
128 
160 

9005 



X 

714 

712 

668 

656 

652 

648 

640 

638 

636 

630 

622 

618 

596 

588 

580 

576 

570 

554 

530 

518 

512 

506 

504 

500 

494 

490 

490 

488 

486 

476 

464 

460 

454 

452 

450 

434 

428 

424 

416 

374 

374 

374 

364 

354 

352 

330 

326 

320 

314 

288 

262 



Y 

521 

357 

415 

320 

477 

519 

351 

503 

454 

292 

392 

310 

361 

434 

344 

408 

311 

S03 

325 

348 

380 

302 

415 

470 

432 

499 

364 

392 

521 

326 

470 

288 

369 

394 

342 

386 

303 

409 

360 

429 

393 

297 

332 

475 

284 

327 

421 

369 

296 

322 

281 



NewX 
696 
694 
650 
662 
634 
630 
646 
620 
638 
636 
624 
624 
602 
590 
586 
578 
576 
556 
536 
524 
518 
512 
506 
502 
496 
492 
496 
490 
488 
482 
466 
466 
456 
454 
456 
436 
434 
426 
418 
376 
376 
376 
366 
356 
354 
332 
322 
316 
316 
290 
264 



NewY 
484 
320 
378 
270 
440 
482 
301 
466 
398 
242 
336 
260 
311 
378 
294 
352 
261 
447 
275 
298 
330 
252 
366 
414 
383 
443 
314 
343 
465 
276 
421 
238 
320 
345 
292 
337 
253 
360 
311 
380 
344 
266 
283 
42 6 
253 
296 
386 
334 
265 
291 
250 



PI 

•1.00 
•1.02 
-1.33 
-1.25 
-1.45 
-1.48 
-1.36 
-1.55 
-1.42 
-1.43 
-1.52 
-1.52 
-1.68 
-1.76 
-1.79 
-1.85 
-1.86 
-2.00 
-2.11 
-2.18 
-2.21 
-2.24 
-2.30 
-2.33 
-2.39 
-2.42 
-2.39 
-2.44 
-2.46 
-2.51 
-2.65 
-2.65 
-2.74 
-2.76 
-2.7 4 
-2.92 
-2.93 
-3.00 
-3.05 
-3.31 
-3.31 
-3.31 
-3.38 
-3.44 
-3.45 
-3.59 
-3.65 
-3 
-3 
-3.84 
-4.00 



.68 
.68 



19.00 

24.15 

21.28 

26.90 

20.48 

19.08 

25.19 

19.68 

21.02 

27.56 

23.27 

27.45 

24.64 

21.28 

25.58 

22.20 

27.40 

20.39 

26.63 

25.36 

23.60 

27.89 

21.43 

20.82 

21.21 

20.45 

24.48 

22.88 

19.72 

26.57 

20.73 

27.34 

24.15 

22.77 

25.69 

23.21 

27.84 

21.50 

24.64 

21.25 

22.83 

27.12 

26.19 

20.66 

27.84 

25.47 

21.17 

23.38 

27.18 

25.75 

28.00 



II 



2.609 
2.703 
0.514 
1.309 
0.173 
0.159 
1.385 
0.105 
0.848 
0.399 
0.282 
0.349 
0.075 
0.426 
0.088 
0.014 
0.007 
2.803 
0.062 
2.230 
0.177 
6.837 
2.542 
0.075 
0.691 
0.031 
1.919 
0.384 
0.278 
1.525 
0.051 
0.083 
1.840 
0.161 
0.484 
0.054 
1.068 
2.862 
0.415 
0.119 
6.077 
0.655 
0.126 
0.250 
0,232 
0.416 
062 
406 
050 
489 
348 



Area 

4.350 
4.640 
2.010 
2.760 
0.870 
0.900 
2.360 
0.580 
2.120 
1.190 
1.350 
1.250 
0.560 
1.830 
0.770 
0.210 
0.130 
3.290 



0 

2, 

0, 
6. 



.290 
.460 
.500 
.150 
3.260 
0.770 
1.720 
0.400 
1.720 
0.690 
1.460 
1.800 
0.690 
0.560 
2.070 
0.370 
1.090 
0.160 
1.750 
3.070 
1.400 
0.720 
5.960 
1.670 
0.820 
1.510 
1.190 
1.400 
0 .740 
1.430 
0.530 
1.560 
4.350 



Ht 

30 
36 
20 
20 
14 
14 
18 
10 
18 
12 
16 
12 

6 
16 

8 

4 

2 
20 

6 
20 

8 
32 
20 
10 
16 

6 
16 

8 
14 
22 

8 

8 
18 

8 
12 

4 
16 
22 
12 

8 
32 

ie 

\2 
U 

12 
14 

10 
14 
8 
16 
26 



Wd 

26 
34 
20 
32 
10 
12 
24 
12 
26 
20 
16 
18 
16 
22 
18 
10 
10 
32 
8 
24 
12 
42 
38 
16 
20 
14 
24 
16 
20 
18 
16 
12 
22 
12 
18 
8 
22 
28 
20 
18 
38 
22 
14 
20 
20 
20 
16 
20 
14 
18 
36 



APPENDIX C5-£ 



WO 91/19274 PCT/US91/03620 

75 

Print; file namg; b.att. Time: 15:14 Date: 12/1/1 , (^0 Page: 

161 262 330 264 299 -4.00 25.30 0.035 0.190 4 8 
9006 244 381 240 346 -4.21 22.71 11.183 9.490 40 54 



APPENDIX C5-3 



WO 91/19274 



76 



PCT/US91/03620 



P^int; fug name: b.ala. T i me: 15:?Q Data- 12/1/10^Q 



Filename: b.als 



Window 


= X : 


244 - 714 








Y : 


281 - 521 






Rec# Spotname_ 


_ Image 


NewX 


NewY 


9001 1 




b 


696 


484 


178 




b 


694 


320 


198 




b 


650 


378 


175 




b 


662 


270 


213 




b 


634 


440 


220 




b 


630 


482 


176 




b 


646 


301 


212 




b 


620 


466 


211 




b 


638 


398 


139 




b 


636 


242 


195 




b 


624 


336 


140 




b 


624 


260 


174 




b 


602 


311 


197 




b 


590 


378 


173 




b 


586 


294 


196 




b 


578 


352 


137 




b 


576 


261 


9002 2 




b 


556 


447 


170 




b 


536 


275 


171 




b 


524 


298 


172 




b 


518 


330 


9003 3 




b 


512 


252 


193 




b 


506 


366 


208 




b 


502 


414 


194 




b 


496 


383 


209 




b 


492 


443 


168 




b 


496 


314 


191 




b 


490 


343 


218 




b 


488 


465 


166 




b 


482 


276 


207 




b 


466 


421 


135 




b 


466 


238 


169 




b 


456 


320 


192 




b 


454 


345 


167 




b 


456 


292 


189 




b 


436 


337 


134 




b 


434 


253 


9004 4 




b 


426 


360 


165 




b 


418 


311 


188 




b 


376 


380 


186 




b 


376 


344 


131 




b 


376 


266 


163 




b 


366 


283 


206 




b 


356 


426 


130 




b 


354 


253 


162 




b 


332 


296 


187 




b 


322 


386 


164 




b 


316 


334 


128 




b 


316 


265 



II 



0. 
0. 
0. 

2, 
0, 



2.609 
2.703 
0.514 
1.309 
0.173 
0.159 
1.385 
0.105 
0.848 
0.399 
0.282 
0.349 
0.075 
0.426 
.088 
,014 
.007 
.803 
.062 
2.230 
0.177 
6.837 
2.542 
0.075 
0.691 
0.031 
1,919 
0.384 
0.278 
1.525 
0.051 
0.083 
1.840 
0.161 
0.484 
0.054 
1.068 
2.862 
0.415 
0.119 
6.077 
0.655 
0.126 
0.250 
0.232 
0.416 
0.062 
0.406 
0.050 



Area 

4.350 
4.640 
2.010 
2.760 
0.870 
0.900 
2.360 
0.580 
2.120 
1.190 
1.350 
1.250 
0.560 
1.830 
0.770 
0.210 
0.130 
3.290 
0.290 
2.460 
0.500 
6.150 
3.260 
0.770 
1.720 
0.400 
1.720 
0.690 
1.460 
1.800 
0.690 
0.560 
2.070 
0.370 
1.090 
0.160 
1.750 
3.070 
1.400 
0.720 
5.96Q 
1.670 
0.820 
1.510 
1.190 
1.400 
0.740 
1.430 
0.530 



Ht 



Wd 



30 
36 
20 
20 
14 
14 
18 
10 
18 
12 
16 
12 

6 
16 

8 

4 

2 
20 

6 
20 

8 
32 
20 
10 
16 

6 
16 

8 
14 
22 

8 

8 
18 

8 
12 

4 
16 
22 
12 

8 
32 
16 
12 
14 
12 
14 
10 
14 

8 



26 
34 
20 
32 
10 
12 
24 
12 
26 
20 
16 
18 
16 
22 
18 
10 
10 
32 

8 
24 
12 
42 
36 
16 
20 
14 
24 
16 
20 
18 
16 
12 
22 
12 
18 

8 
22 
28 
20 
18 
38 
22 
14 
20 
20 
20 
16 
20 
14 



APPENDIX C5-<* 



WO 91/19274 PCIYUS91/03620 

77' 

Print file name: b.ala. Time: 15:20 Date: 1 2/1/1989 Page: 



160 




b 


290 


291 


0-489 


1.560 


16 


18 


9005 


5 


b 


264 


250 


4.348 


4.350 


26 


36 


161 




b 


264 


299 


0.035 


0.190 


4 


8 


9006 


6 


b 


240 


346 


11.183 


9.490 


40 


54 



APPENDIX C5-5 



WO 91/19274 PCT/US91/03620 

78 

Print fU« aamSJ b,mks. Time: 15:19 p a i- 0 . i^y io.n PaQe ■ 

b.mks 



Window = x : 


244 


- 714 


Y : 


281 


~ 521 




6 Markers 


Ma rke rename 


X 


y 


9001 


714 


521 


9002 


554 


503 


9003 


506 


302 


9004 


424 


409 


9005 


262 


281 


9006 


244 


381 



APPENDIX C5-6 



WO 91/19274 



79 



PCT/US91/03620 




APPENDIX C6-1 



WO 91/19274 



do 



Print £41s name: c.att. Time: 15:15 n ate; 12/yiQRQ 
Filename: c.att 



PCT/US91/03620 



Pace : 



Rec# 

9001 
120 
168 
177 
166 
148 
167 
176 
117 
175 
145 
165 
146 
116 
164 

9002 
141 
162 

9003 
171 
163 
140 
173 
142 
144 
114 
172 
143 
138 
112 

9004 
139 
136 
158 
109 
160 
135 
134 
159 
137 
106 
104 
133 
132 

9005 

9006 
131 



682 
660 
678 
674 
674 
672 
648 
640 
636 
634 
622 
622 
606 
600 
582 
554 
508 
502 
500 
498 
488 
486 
484 
484 
480 
470 
464 
448 
446 
426 
426 
416 
376 
376 
374 
374 
372 
330 
328 
320 
318 
316 
288 
260 
258 
242 
242 



514 
295 
412 
475 
386 
352 
403 
470 
310 
508 
346 
440 
382 
301 
423 
486 
337 
402 
294 
455 
421 
333 
508 
353 
379 
316 
457 
363 
336 
301 
404 
358 
352 
393 
300 
428 
347 
335 
427 
377 
308 
303 
336 
341 
297 
397 
361 



NewX 
696 
692 
692 
688 
688 
686 
662 
654 
648 
648 
634 
624 
608 
612 
584 
556 
520 
502 
512 
500 
488 
498 
486 
496 
480 
482 
4 64 
448 
458 
438 
426 
416 
376 
376 
380 
374 
372 
336 
326 
318 
324 
322 
294 
266 
264 
240 
240 



NewY 
484 
253 
382 
445 
356 
322 
373 
440 
268 
478 
304 
401 
343 
259 
384 
447 
295 
358 
252 
416 
377 
291 
469 
311 
335 
274 
413 
319 
294 
259 
360 
314 
308 
349 
253 
384 
303 
268 
376 
326 
261 
256 
289 
294 
250 
346 
310 



— i 

-1 

-1 

•1 

-1 

-1 

-1 

-1 

-1 

-1 

-1 

-1 

-1 

-1 

-1 

-2 

-2 

-2 

-2 

-2 

-2 

-2 

-2 

-2 

-2 

-2 

-2 

-2 

-2 

-2 

-3 

-3 

-3 

-3 

-3 

-3 

-3 

-3 

-3 

-3 

-3 

-3 

-3 

-3 

-4 

-4 

-4 



MW 

00 

03 

03 

06 

06 

08 

25 

30 

35 

35 

45 

52 

63 

60 

80 

00 

20 

33 

24 
,35 
,46 
,37 

47 
,39 
.53 
,51 
,67 
,81 
.72 
,90 
.00 
.07 
.31 
.31 
,29 
.33 

34 
.56 



II 



.62 
• 67 
.63 
.65 
.82 
.99 
.00 
.21 
.21 



19.00 
27.84 
21.22 
20.42 
21.85 
24.04 
21.34 
20.48 
27.01 
19.23 
25.03 
20.98 
22.88 
27.51 
21.20 
20.39 
25.53 
21.68 
27.89 
20.79 
21.29 
25.75 
19.57 
24.64 
23.32 
26.68 
20.83 
24.20 
25.58 
27.51 
21.50 
24.48 
24.81 
22.46 
27.84 
21.20 
25.08 
25.91 
21.30 
23.82 
27.40 
27.67 
25.86 
25.58 
28.00 
22.71 
24.70 



1. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 

1. 

0. 
0. 
0. 
0. 

1. 
1. 

2. 
5. 
0. 
0. 
0. 
0. 
0. 
0. 

1. 

0. 
0. 
0. 
0. 
4. 
0. 
0. 
5. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0, 
2, 
9. 
0, 



216 
012 
131 
148 
076 
360 
286 
289 
713 
062 
103 
462 
114 
425 
402 
703 
218 
658 
686 
029 
904 
172 
413 
659 
165 
003 
067 
519 
427 
747 
198 
485 
026 
207 
403 
121 
012 
608 
060 
582 
019 
006 
707 
013 
975 
890 
486 



Area 

2.860 
0.080 
0.720 
1.010 
0.480 
1.620 
1.170 
1.270 
1.830 
0.420 
2.040 
1.800 
0.850 
1.430 
1.880 
2.620 
1.860 
3.580 
5.700 
0.370 
1.880 
0.290 
1.910 
0.980 
0.240 
1.350 
0.660 
1.220 
1.110 
1.780 
4.060 
1.430 
0.160 
5.860 
1.540 
0.770 
0.110 
1.800 
0.720 
1.700 
0.240 
0.080 
1.930 
0.080 
3.900 
9.940 
1.110 



Ht 



Wd 



20 

2 
12 
12 

8 
16 
14 
14 
16 

8 
20 
18 
12 
14 
18 
20 
18 
22 
34 

6 
18 

6 
18 
10 

6 
18 
10 
14 
12 
18 
30 
16 

4 
30 
16 

8 

2 
18 
10 
18 

4 

2 
20 

? 
26 
36 
10 



26 

6 
16 
18 
12 
18 
18 
18 
26 
10 
20 
20 
12 
20 
20 
26 
20 
34 
36 
12 
20 

8 
22 
18 
10 
16 
14 
14 
16 
18 
28 
16 

8 
36 
18 
18 

8 
18 
14 
20 
10 

6 
18 

6 
34 
54 
22 



APPENDIX C6-2 



WO 91/19274 

Si 

Print; file nam$; c T al g> Time ; 1?;2Q Date: 12/1/19?*° 



PCT/US91/03620 

Page; . 



Filename: c.als 
Window = X : 242 - 682 
Y : 294 ~ 514 



Rec# 


Spot name 


Image 


NewX 


A A f\ 1 

9001 


l 


c 


r- 1\ r 

696 


120 




c 


692 


168 




c 


692 


177 




c 


688 


166 




c 


688 


148 




c 


686 


167 




c 


662 


176 




c 


654 


117 




c 


648 


175 




c 


648 


145 




c 


634 


165 




c 


624 


146 




c 


608 


116 




c 


612 


164 




c 


584 


9002 


2 


c 


556 


141 




c 


520 


162 




c 


502 


9003 


3 


c 


512 


171 




c 


500 


163 




c 


488 


140 




c 


498 


173 




c 


466 


142 




c 


496 


144 




c 


480 


114 




c 


482 


172 




c 


464 


143 




c 


448 


138 




c 


458 


112 




c 


438 


9004 


4 


c 


426 


139 




c 


416 


136 




c 


376 


158 




c 


376 


109 




c 


380 


160 




c 


374 


135 




c 


372 


134 




c 


336 


159 




c 


326 


137 




c 


318 


106 




c 


324 


104 




c 


322 


133 




c 


294 


132 




c 


266 


9005 


5 


c 


264 


9006 


6 




240 


131 




c 


240 



NewY II Area Ht wd 



484 


1.216 


2 . 860 


20 


26 


253 


0 . 012 


0 . 080 


2 


6 


382 


0.131 


0.720 


12 


16 


445 


0.148 


1.010 


12 


18 


356 


0.076 


0.480 


8 


12 


322 


0.360 


1.620 


16 


18 


373 


0 .286 


1.170 


14 


18 


440 


0.289 


1.270 


14 


18 


268 


0.713 


1.830 


16 


26 


478 


0.062 


0.420 


8 


10 


304 


1.103 


2.040 


20 


20 


401 


0.462 


1.800 


18 


20 


343 


0.114 


0.850 


12 


12 


259 


0.425 


1.430 


14 


20 


384 


0.402 


1.880 


18 


20 


447 


1.703 


2.620 


20 


26 


295 


1.218 


1.860 


18 


20 


358 


2. 658 


3.580 


22 


34 


252 


5.686 


5.700 


34 


36 


416 


0.029 


0.370 


6 


12 


377 


0.904 


1.880 


18 


20 


291 


0.172 


0.290 


6 


8 


469 


0.413 


1.910 


18 


22 


311 


0.659 


0.980 


10 


18 


335 


0.165 


0.240 


6 


10 


274 


1.003 


1.350 


18 


16 


413 


0.067 


0.660 


10 


14 


319 


0.519 


1.220 


14 


14 


294 


0.427 


1.110 


12 


16 


259 


0.747 


1.780 


18 


18 


360 


4.198 


4.060 


30 


28 


314 


0.485 


1.430 


16 


16 


308 


0.026 


0.160 


4 


8 


349 


5.207 


5.860 


30 


36 


253 


0.403 


1.540 


16 


18 


384 


0.121 


0.770 


8 


18 


303 


0.012 


0.110 


2 


8 


288 


0.608 


1.800 


18 


18 


376 


0.060 


0.720 


10 




326 


0.582 


1.700 


18 


20 


261 


0.019 


0.240 


4 


10 


256 


0.006 


0.080 


a 


6 


289 


0.707 


1.930 


20 


18 


294 


0,013 


o.oao 


a 


6 


250 


2.975 


3.900 


26 


34 


346 


9.890 


9.940 


36 


54 


310 


0.466 


1.1x0 


10 


22 



APPENDIX C6-3 



WO 91/19274 



PCT/US91/03620 



82 

Print file name; c.mks. Time: 15-1 9 Date; 12/1/lQflQ 
c .xiUcs 

Window = X : 242 - 682 
Y : 294 ~ 514 



Marker name 



X 



6 Markers 
Y 



9001 
9002 
9003 
9004 
9005 
9006 



682 514 

554 486 

500 294 

426 404 

258 297 

242 397 



APPENDIX 



WO 91/19274 PCT/US91/03620 

Print; fi,le aaffiSJ abc.iumh, Time: 17:59 bate : 3/8/1990 Page: 

#NS 3 /* # of scans in the match group */ 
#PI 0.16 /* PI maximum tolerance */ 
#MW 1.12 /* MW maximum tolerance */ 
I;C:I:I:R:R:R:R:I:I:I:I 



REC 


GEL 


X 


Y 


PI 


MW 




II 


AREA 


MAT a 




b 


9001 


a 


696 


484 


-1.00 


19.00 


3 


.652 


5.140 


o 

■J 




9001 


9001 


218 


b 


486 


521 


-2.46 


19.71 


G 


.278 


1.460 




ICQ 


218 


173 


166 


b 


476 


326 


-2.50 


26.56 


1 


.525 


1.800 






166 


114 


113 


a 


690 


319 


-1.04 


24.20 


3 


.868 


5.060 


"3 


±± o 


178 


148 


137 


a 


658 


369 


-1.27 


21.38 


0 


.859 


1.990 


J 


x O l 


198 


167 


112 


a 


650 


274 


-1.33 


26. 68 




.391 


1.670 


•a 


1 1 o 
lie 


175 


117 


157 


a 


636 


457 


-1.42 


20.02 


1 


.162 


2.200 


■a 
J 


1 C"7 


212 


175 


110 


a 


636 


308 


-1.42 


24.80 


2 


.391 


3.340 


j 


11U 


176 


145 


149 


a 


630 


402 


-1.48 


20.96 


1 


.845 


2.860 


•a 


1 A Q 


211 


165 


207 


b 


464 


470 


-2,65 


20.72 


0 


.051 


0.690 




145 


207 


172 


135 


a 


620 


348 


-1.54 


22.54 


0 


.281 


1.140 


o 


135 


195 


146 


107 


a 


620 


268 


-1.54 


27.01 


0 


.049 


0.190 


3 


107 


140 


116 


167 


b 


450 


342 


-2.74 


25.69 


0 


.484 


1.090 


3 


103 


167 


138 


148 


a 


586 


389 


-1.78 


21.12 


0 


.264 


1.540 


3 


148 


197 


164 


9006 


c 


242 


397 


-4.21 


22.70 


9 


.890 


9.940 


3 


9006 


9006 


9006 


159 


c 


328 


427 


-3.61 


21.29 


0 


.060 


0.720 


3 


12 9 


187 


159 


9005 


c 


258 


297 


-4.00 


28.00 


2 


.975 


3.900 


3 


ft A ft C 

9005 


9005 


9005 


9002 


a 


556 


447 


-2.00 


20.38 


4 


.757 


5.120 


3 


9002 


9002 


9002 


106 


a 


524 


302 


-2.18 


25.13 


1.773 


2.330 


3 


106 


171 


141 


9003 


a 


512 


252 


-2.24 


27.88 


8 


.928 


6.710 


3 


9003 


9003 


9003 


133 


a 


504 


372 


-2.31 


21.35 


3 


.676 


5.800 


3 


133 


193 


162 


146 


a 


500 


417 


-2.34 


20.78 


0 


.203 


1.300 


3 


146 


208 


171 


155 


a 


492 


464 


-2.42 


19.76 


0 


.379 


1.910 


3 


155 


218 


173 


104 


a 


490 


319 


-2.44 


24.20 


1 


.404 


1.380 


3 


104 


168 


142 


132 


a 


480 


333 


-2.52 


23.43 


0 


.039 


0.080 


3 


132 


191 


144 


102 


a 


478 


277 


-2.54 


26.52 


1.250 


1.540 


3 


102 


166 


114 


145 


a 


460 


417 


-2.70 


20.78 


0 


.140 


1.140 


3 


145 


207 


172 


103 


a 


454 


298 


-2.75 


25.36 


0 


.437 


0.980 


3 


103 


167 


138 


131 


a 


452 


324 


-2.77 


23.93 


1 


.888 


2.170 


3 


131 


169 


143 


100 


a 


432 


261 
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22.45 
27.51 
27.84 
27.12 
24.30 
24.47 
24.63 
21.50 



II 

0.148 
1.023 
0.399 
1.006 
0.075 
0.022 
0.088 
0.111 
0.691 
0.904 
0.035 
0.013 
0.076 
0.300 
0.289 
0.173 
0.052 
0.007 
0.531 
0.486 
0.296 
0.707 
0.489 
0.150 
0.050 
0.019 
0.191 
0.406 
0.582 
0.162 
0.416 
0.608 
0.026 
0.012 
0.126 
0.051 
0.121 
0.119 
3.781 
6.077 
5.207 
0 .702 



403 
655 
131 
485 
415 



4.162 



AREA 

1.010 
2.620 
1.190 
1.800 
0.560 
0.190 
0.770 
1.090 
1.720 
1.880 
0.190 
0.080 
0.480 
1.030 
1.270 
0.870 
0.720 
0.130 
1.190 
1.110 
1.540 
1.930 
1.560 
0.980 
0.530 
0.240 
1.300 
1.430 
1.700 
1.110 
1.400 
1.800 
0.270 
0.110 
0.820 
0.370 
0.770 
0.720 
5.190 
5.960 
5.860 
1.960 
1.540 
1.670 
0.610 
1.430 
1.400 
3.900 



MAT a 



150 
150 
73 
73 
109 
109 
108 
108 
0 
0 
0 
0 

136 
136 
0 
0 

105 
105 
94 
94 
95 
95 
95 
96 
96 
96 
127 
127 
3 127 
3 99 
3 99 
3 99 
3 98 
3 98 
3 98 
3 144 
3 144 
3 144 
3 128 
3 128 
3 128 
3 97 
3 97 
3 97 
3 101 
3 101 
3 101 
3 9004 



0 
0 

139 
139 
174 
174 
173 
173 
194 
194 
161 
161 
0 
0 

213 
213 
137 
137 
0 
0 

160 
160 
160 
128 
128 
128 
164 
164 
164 
162 
162 
162 
163 
163 
163 
188 
188 
188 
186 
186 
186 
131 
131 
131 
165 
165 
165 
9004 



177 
177 
0 
0 
0 
0 
0 
0 

163 
163 
132 
132 
166 
166 
176 
176 
0 
0 

131 
131 
133 
133 
133 
106 
106 
106 
137 
137 
137 
134 
134 
134 
135 
135 
135 
160 
160 
160 
158 
158 
158 
109 
109 
109 
139 
139 
139 
9004 
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9004 


b 


424 


409 


-3 . 00 


21.50 


2.862 


3.070 


3 


9004 


9004 


9004 


9004 


c 


426 


404 


-3.00 


21.50 


4.198 


4.060 


3 


9004 


9004 


9004 


100 


a 


432 


261 


-2.95 


27.39 


0.884 


1.910 


3 


100 


134 


112 


112 


c 


426 


301 


-2.90 


27.51 


0.747 


1.780 


3 


100 


134 


112 


134 


b 


428 


303 


-2.93 


27.84 


1.068 


1.750 


3 


100 


134 


112 


131 


a 


452 


324 


-2.77 


23.93 


1.888 


2.170 


3 


131 


169 


143 


143 


c 


448 


363 


-2.80 


24.20 


0.519 


1.220 


3 


131 


169 


143 


169 


b 


454 


369 


-2.74 


24.14 


1.840 


2.070 


3 


131 


*69 


143 


132 


a 


480 


333 


-2.52 


23.43 


0.039 


0.080 


3 


132 


191 


144 


191 


b 


488 


392 


-2.44 


22.87 


0.384 


0.690 


3 


132 


191 


144 


144 


c 


480 


379 


-2.52 


23.31 


0.165 


0.240 


3 


132 


191 


144 


104 


a 


490 


319 


-2.44 


24.20 


1 . 404 


1.380 


3 


104 


168 


142 


168 


b 


490 


364 


-2.39 


24.47 


1.919 


1.720 


3 


104 


168 


142 


142 


c 


484 


353 


-2.39 


24.63 


0.659 


0.980 


3 


104 


168 


142 


146 


a 


500 


417 


•2.34 


20.78 


0.203 


1.300 


3 


146 


208 


171 


208 


b 


500 


470 


-2.32 


20.81 


0.075 


0.770 


3 


146 


208 


171 


171 


c 


498 


455 


-2.34 


20.79 


0.029 


0.370 


3 


146 


208 


171 


133 


a 


504 


372 


-2.31 


21.35 


3.676 


5.800 


3 


133 


193 


162 


193 


b 


504 


415 


-2.29 


21.43 


2.542 


3.260 


3 


133 


193 


162 


162 


c 


502 


402 


-2.32 


21.68 


2.658 


3.580 


3 


133 


193 


162 


9003 


a 


512 


252 


-2.24 


27.88 


8.928 


6.710 


3 


9003 


9003 


9003 


9003 


b 


506 


302 


-2.24 


27.88 


6.837 


6.150 


3 


9003 


9003 


9003 


9003 


c 


500 


294 


-2.24 


27.88 


5.686 


5.700 


3 


9003 


9003 


9003 


106 


a 


524 


302 


-2.18 


25.13 


1.773- 


2.330 


3 


106 


171 


141 


171 


b 


518 


348 


-2.18 


25.36 


2.230 


2.460 


3 


106 


171 


141 


141 


c 


508 


337 


-2.20 


25.53 


1.218 


1.860 


3 


106 


171 


141 


9002 


a 


556 


447 


-2.00 


20.38 


4.757 


5.120 


3 


9002 


9002 


9002 


9002 


b 


554 


503 


-2.00 


20.38 


2.603 


3.290 


3 


9002 


9002 


9002 


9002 


c 


554 


486 


-2.00 


20.38 


1.703 


2.620 


3 


9002 


9002 


9002 


9005 


c 


258 


297 


-4.00 


28.00 


2.975 


3.900 


3 


9005 


9005 


9005 


9005 


a 


264 


250 


-4.00 


28.00 


3.127 


4.590 


3 


9005 


9005 


9005 


9005 


b 


262 


281 


-4.00 


28.00 


4.348 


4.350 


3 


9005 


9005 


9005 


159 


c 


328 


427 


-3.61 


21.29 


0.060 


0.720 


3 


129 


187 


159 


129 


a 


332 


382 


-3.58 


21.21 


0.074 


0.930 


3 


129 


187 


159 


187 


b 


326 


421 


-3.65 


21.17 


0.062 


0.740 


3 


129 


187 


159 


9006 


c 


242 


397 


-4.21 


22.70 


9.890 


9.940 


3 


9006 


9006 


9006 


9006 


a 


240 


346 


-4.21 


22.70 


7.771 


7.790 


3 


9006 


9006 


9006 


9006 


b 


244 


381 


-4.21 


22.70 


11.183 


9.490 


3 


9006 


9006 


9006 


148 


a 


586 


369 


-1.78 


21.12 


0.264 


1.540 


3 


148 


197 


164 


197 


b 


588 


434 


-1.75 


21.28 


0.426 


1.830 


3 


148 


197 


164 


164 


c 


582 


423 


-1.79 


21.20 


0.402 


1.880 


3 


148 


197 


164 


167 


b 


450 


342 


-2.74 


25.69 


0.484 


1.090 


3 


103 


167 


138 


103 


a 


454 


298 


-2.75 


25.36 


0.437 


0.980 


3 


103 


167 


138 


138 


c 


446 


336 


-2.72 


25.57 


0.427 


1.110 


3 


103 


167 


138 


107 


a 


620 


268 


-1.54 


27.01 


0.049 


0.190 


3 


107 


140 


116 


140 


b 


618 


310 


-1.51 


27 .45 


0.349 


1.250 


3 


107 


140 


116 


116 


c 


600 


301 


-1.60 


27.51 


0.425 


1.430 


3 


107 


140 


116 


135 


a 


620 


348 


-1.54 


22.54 


0.281 


1.140 


3 


135 


195 


146 


195 


b 


622 


392 


-1.51 


23.27 


0.282 


1.350 


3 


135 


195 


146 


146 


c 


606 


382 


-1.62 


22.87 


0.114 


0.850 


3 


135 


195 


146 


207 


b 


464 


470 


-2.65 


20.72 


0.051 


0.690 


3 


145 


207 


172 


145 


a 


460 


417 


-2.70 


20.78 


0.140 


1.140 


3 


145 


207 


172 


172 


c 


464 


457 


-2.67 


20.82 


0.067 


0.660 


3 


145 


207 


172 


149 


a 


630 


402 


-1.48 


20.96 


1.845 


2.860 


3 


149 


211 


165 
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Time; 11:06 Dajgj 1 2/5/lQfiQ 



211 b 

165 c 
110 a 
176 
145 
112 
175 
117 
137 
198 
167 c 

113 a 
178 b 
148 c 

166 b 
102 a 

114 c 
218 b 
155 a 
173 

9001 
9001 
9001 



636 


454 


-1. 


41 


21 


.02 


0.848 


2 120 




1 A Q 


Oil 

211 


165 


622 


440 


-1. 


51 


20 


.97 


0 . 462 


1.800 


o 






165 


636 


308 


-1. 


42 


24 


. 80 


2 3Q1 






110 


176 


145 


640 


351 


-1. 


36 


25 


. 19 


" 3R ^ 


2 *3 c n 


3 


110 


176 


145 


622 


346 


-1. 


45 


25 


. 03 




2 . U4U 


3 
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145 
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-1. 


33 


26 
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i 3 q n 
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3 
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35 


27 


01 
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27 


2 1 
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1 . 990 


3 


137 


198 


167 
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415 


-1. 


33 


21 


2 R 




2 .010 


3 


137 


198 


167 


646 


403 


-1.25 


21 


34 


n 9QC 
U . Z O D 


1 . 170 


3 


137 


198 


167 


690 


319 


-1. 


04 


24 


2 n 


O.ODO 


5 .060 


3 


113 


178 


148 


712 


357 


-1. 


01 


24 


1 4 




4 . 640 


3 


113 


178 


148 


672 


352 


-1. 


08 


24 


.04 




1 . 620 


3 


113 


178 


148 


476 


326 


-2. 


50 


26 


.56 


1 .525 


•i • S> \J V 






loo 


114 


478 


277 


-2. 


54 


26 


.52 


1.250 


1.540 


3 


102 


166 


114 


470 


316 


-2. 


50 


26 


.68 


1.003 


1.350 


3 


102 


166 


114 


486 


521 


-2. 


46 


19 


.71 


0.278 


1.460 


3 


155 


218 


173 


492 


464 


-2. 


42 


19 


.76 


0.379 


1.910 


3 


155 


218 


173 


484 


508 


-2. 


47 


19 


.56 


0.413 


1.910 


3 


155 


218 


173 


696 


484 


-1. 


00 


19 


.00 


3.652 


5.140 


3 


9001 


9001 


9001 


714 


521 


-1. 


00 


19 


.00 


2.609 


4.350 


3 


9001 


9001 


9001 


682 


514 


-1. 


00 


19 


.00 


1.216 


2.860 


3 


9001 


9001 


9001 
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#NS 3 /* # of scans in the match group */ 
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IEC 


GEL 


X 


Y 


PI 


MW 


II 


AREA 


MAT a 




b 


-10 


abc 


-1 


-1 


-4.21 


24.65 


0.508 


1.150 


Cm 


Q a 


0 


131 


-6 


abc 


-1 


-1 


-3.99 


25.43 


0.024 


0.135 


o 




161 


132 


-11 


abc 


-1 


-1 


-3.82 


25.80 


0.497 


1.677 


5 


95 


160 


133 


-13 


abc 


-1 


-1 


-3.66 


23.55 


0.393 


1.477 


3 


127 


164 


137 


104 


c 


316 


303 


-3.65 


27.67 


0.006 


0.080 


1 


0 


0 


104 


-12 


abc 


-1 


-1 


-3.65 


27.32 


0.073 


0.583 


3 


96 


128 


106 


-31 


abc 


-1 


-1 


-3.61 


21.22 


0.065 


0.797 


3 


12 9 


187 


159 


-14 


abc 


-1 


-1 


-3.56 


25.70 


0.395 


1.437 


3 


99 


162 


134 


130 


b 


352 


284 


-3.45 


27.84 


0.232 


1.190 


1 


0 


130 


0 


206 


b 


354 


475 


-3.44 


20.65 


0.250 


1.510 


1 


0 


206 


0 


-15 


abc 


-1 


-1 


-3.36 


25.61 


0.055 


0.400 


3 


98 


163 


135 


-16 


abc 


-1 


-1 


-3.32 


21.21 


0.097 


0.620 


3 


144 


188 


160 


136 


c 


376 


352 


-3.30 


24.80 


0.026 


0.160 


1 


0 


0 


136 


-17 


abc 


-1 


-1 


-3.30 


22.57 


5.022 


5.670 


3 


128 


186 


158 


-18 
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-1 


-1 


-3.29 


27.49 


0.587 


1.723 


3 


97 


131 


109 


-19 


abc 


-1 


-1 


-3*05 


24.47 


0.344 


1.147 


3 
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165 
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-21 
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-1 


-1 


-2.93 


27.58 


0.900 


1.813 


3 
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134 
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189 


b 


434 


386 


-2.92 


23.20 


0.054 


0.160 


1 


0 
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0 


-22 


abc 


-1 


-1 


-2.77 


24.09 


1.416 


1.820 


3 


131 


169 


143 


192 


b 


452 


394 


-2.75 


22.77 


0.161 


0.370 


1 


0 


192 


0 


-34 


abc 


-1 


-1 


-2.74 


25.54 


0.449 


1.060 


3 


103 


167 


138 


-37 


abc 


-1 


-1 


-2.67 


20.77 


0 .086 


0 .830 


3 


145 


207 


172 


135 


b 


460 


288 


-2.65 


27.34 


0.083 


0.560 


1 


0 


135 


0 


-45 


abc 


-1 


-1 


-2.51 


26.59 


1.259 


1.563 


3 


102 


166 


114 


-23 


abc 


-1 


-1 


-2.49 


23.20 


0 . 196 


0.337 


3 


132 


191 


144 


-46 


abc 


-1 


-1 


-2.45 


19.68 


0.357 


1.760 


3 


155 


218 
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-5 
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-1 


-1 


-2.43 


21.25 


0.798 


1.800 


2 


0 


194 


163 


209 


b 


490 


499 


-2.42 


20.45 


0.031 


0.400 


1 


0 


209 


0 


-24 


abc 


-1 


-1 


-2.41 


24.43 


1.327 


1.360 


3 


104 


168 


142 


140 


c 


486 


333 


-2.36 


25.75 


0.172 


0.290 


1 


0 


0 


140 


-25 


abc 


-1 


-1 


-2.33 


20.79 


0.102 


0.813 


3 


14 6 


208 


171 


-26 


abc 


-1 


-1 


-2.31 


21.49 


2.959 


4.213 


3 


133 


193 


162 


172 


b 


512 


380 


-2.21 


23.60 


0.177 


0.500 


1 


0 
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0 


-28 


abc 


-1 


-1 


-2.19 


25.34 


1.740 


2.217 


3 
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171 


141 


170 


b 


530 


325 


-2.10 


26.62 


0.062 


0.290 


1 


0 


170 


0 


134 


a 


560 


336 


-1.98 


23.27 


0.055 


0.480 


1 


134 


0 


0 


-9 


abc 


-1 


-1 


-1.88 


27.28 


0.030 


0.425 


2 


105 


137 


0 


196 


b 


576 


408 


-1.85 


22.20 


0.014 


0.210 


1 


0 


196 


0 


-4 


abc 


-1 


-1 


-1.82 


25.35 


0.100 


0.930 


2 


108 


173 


0 


-33 


abc 


-1 


-1 


-1.77 


21.20 


0.364 


1.750 


3 


148 


197 


164 


-3 


abc 


-1 


-1 


-1.68 


24.50 


0.049 


0.375 


2 


109 


174 


0 


-36 


abc 


-1 


-1 


-1.56 


22.89 


0.226 


1.113 


3 


135 


195 


146 


-35 


abc 


-1 


-1 


-1.55 


27.32 


0.274 


0.957 


3 


107 


140 


116 


-40 


abc 


0 


0 


-1.48 


19.85 


0.633 


1.390 


2 


157 


212 


0 


-38 


abc 


-1 


-1 


-1.47 


20.98 


1.052 


2.260 


3 


149 


211 


165 


-2 


abc 


-1 


-1 


-1.45 


27.77 


0 .702 


1.495 


2 


73 


139 


0 


-39 


abc 


-1 


-1 


-1.41 


25.01 


1.626 


2.580 


3 
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176 


145 


-41 


abc 


0 


0 


-1.41 


19.15 


0.110 


0.660 
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0 
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Print file nam^ : abe.eim. 
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.Time; 11:05 Da**- 12/S/iQftQ 



-8 abc 
-42 abc 
-43 abc 
-44 abc 

-1 abc 



168 
120 



c 
c 

-7 abc 
-32 abc 
-30 abc 
-20 abc 
-27 abc 
-29 abc 
-47 abc 



-1 
-1 
-1 
-1 
-1 
678 
680 
-1 
-1 
-1 
-1 
-1 
-1 
-1 



-1 
-1 
-1 
-1 
-1 
412 
295 
-1 
-1 
-1 
-1 
-1 
-1 
•1 



-1.37 
-1.31 
-1.28 
-1.04 
-1.03 
-1.02 
-1.02 
-1.02 
-4.21 
-4.00 
-3.00 
-2.24 
-2.00 
-1.00 



20.47 
26.86 
21.33 
24.13 
20.51 
21.21 
27.84 
21.81 
22.70 
28.00 
21.50 
27.88 
20.38 
19.00 



0.231 
1.138 
0.553 
2.310 
0.586 
0.131 
0.012 
0.188 
9.615 
3.483 
3.741 
7.150 
3.088 
2.492 



1.070 
2.087 
1.723 
3.773 
1.815 
0.720 
0.080 
0.755 
9.073 
4.280 
3.677 
6.187 
3.677 
4.117 



2 
3 
3 
3 
2 
1 
1 
2 



0 

112 
137 
113 
150 
0 
0 

136 
3 9006 
3 9005 
3 9004 
3 9003 
3 9002 
3 9001 



213 
175 
198 
178 
0 
0 
0 
0 

9006 
9005 
9004 
9003 
9002 
9001 



176 
117 
167 
148 
177 
168 
120 
166 
9006 
9005 
9004 
9003 
9002 
9001 
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Time: 17:33 'Data: 3/6/1990 



PCT/US91/03620 

Page: 1 
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We Claim: 

SO 1. A method of matching two-dimensional (2-D> patterns, 
said method comprising the steps of: 

(a) image scanning a plurality oT 2— D patterns and 
producing corresponding data -files -for each scanned 2— D 
pattern and recording logistic data and pattern physical 

25 characteristics data in a scan coordinate system; 

(b > identifying pattern members in each of said 
scanned 2— D pattern that bound investigative patterns and 
designating said identified pattern members as marker 
members ; 

30 (c> designating one of said plurality of 2-D pattern 

as a reference 2-D pattern, said reference 2-D pattern 
having a reference pattern member data file which includes a 
set of reference marker members and unknown reference 
pattern members each having respective coordinates in said 

35 scan coordinate system; 

(d) designating at least one of remaining ones of said 
plurality of 2-D patterns as a study 2-D pattern, said study 
2-D pattern having a study pattern data file which includes 
a set of study marker members and unknown study pattern 

40 members each having respective coordinates in said scan 
coordinate system; 

(e) performing a first transformation step that 
transforms positional coordinates of said set of reference 
marker members, said unknown reference pattern members and 

45 each member of said set of study marker pattern members, 
from said scan coordinate system to a reference coordinate 
system, said first transformation step resulting in each 
member of said set of study marker members being in a 
registered relationship with a corresponding member of said 

50 set of reference marker members; 

Cf ) performing a second transformation step that 

transforms positional coordinates of each of said unknown 
study pattern members from said scan coordinate system to 
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said reference coordinate system and generating at least one 
adjusted data -file containing new coordinate information -for 
all pattern members manipulated by said first and second 
transformation steps; 
5 (g> repeating said steps <d), (e) and <f> on any 

remaining one of said plurality of S-D pattern; and 

(h) comparing coordinates of all pattern members in 
said at least one adjusted data file and producing results 
identifying matching pattern members. 

10 A method of matching two-dimensional (2-D) patterns, as 

recited in claim 1, wherein said step (h) further includes 
the steps of: 

verifying that potentially matching unknown spots are 
within an acceptable vector area formed by marker spot 

15 vectors before producing said matching results; and 

resolving matching results having contradictory 
matching information among spot pattern members* 
3. A method of matching two-dimensional <E-D> patterns, as 
recited in claim 1, wherein: 

SO said a— D patterns being protein spot patterns in a 

plurality of two— d i mens ional gel electrophoretograms <E-D 
gels) and said step of identifying said pattern members as 
marker members further includes assigning isoelectric 
focusing (PI), and molecular weight CMW) dimensional 

£5 separation values to said marker members; and 

said method further includes after performing said 
second transformation step, interpolating said assigned PI 
and MW separation values to said marker members to determine 
PI and MW values for said unknown study patterns. 
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A method of matching two-dimensional (2-D) patterns, as 
recited in claim 3, wherein: 

said step of comparing includes the step of verifying 
that potentially matching unknown spots are within an 
acceptable vector area formed by marker spot vectors before 
producing said matching results; 

said method further includes the step of analyzing said 
matching results using PI and MW separation values resulting 
from said interpolating step; and 

resolving matching results having contradictory 
matching information among pattern members. 



5- A method of matching two-dimensional CE-D) patterns, as 
recited in claim 3, wherein after said step of comparing and 
producing matching results includes: 

analyzing said matching results using PI and MW 
separation values resulting from said interpolating step; 
and 

repeating said steps (b) through (h) and said analyzing 

step . 



in 



20 6 - A method of matching protein spot patterns 

plurality of two-dimensional gel electrophoretograms <E-D 

gels) » said method comprising the steps of: 

<a) image scanning each of said plurality of E-D gels 

and producing corresponding spot data files for each S-D gel 
25 that contain logistic data and spot physical characteristics 

data in a scan coordinate system; 

(b) identifying spot members in each of said plurality 
of E-D gels that bound investigative spot patterns and 
designating said identified spot members as marker spot 

30 members; 

(c) designating one of said plurality of E-D gels as a 
reference E-D gel , said reference E-D gel having an 
associated reference spot data file including a set of 
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reference marker spot members and unknown reference spot 
members having respective coordinates in said scan 
coord i nate system ; 

<d> designating at least one of remaining ones of said 
5 plurality of E-D gels as a study 2-D gel, said study 2-D gel 
having an associated study spot data file including a set of 
study marker spot members and unknown study spot members 
having respective coordinates in said scan coordinate 
system; 

10 (e) performing a first transformation step that 

transforms positional coordinates of said set of reference 
marker spot members , said unknown reference spot members and 
each member of said set of study marker spot members from 
said scan coordinate system to a reference coordinate 

15 system, said first transformation step resulting in each 
member of said set of study marker spot members being in a 
registered relationship with a corresponding member of said 
set of reference marker spot members; 

<f > performing a second transformation step that 

20 transforms positional coordinates of each of said unknown 
study spot members from said scan coordinate system to said 
reference coordinate system and generating at least one 
adjusted data file containing new coordinate information for 
all spot pattern members manipulated by said first and 

25 second transformation steps; 

(g) repeating said steps (d>, <e> and (f) on any 
remaining one of said plurality of 2— D gels; and 

(h) comparing coordinates of all spot pattern members 
in said at least one adjusted data file and producing 

30 results identifying matching spot pattern members. 

7. A method of matching protein spot patterns in a plurality 
of two— dimensional gel elec trophoretograms 9 as recited in 
claim 6, wherein said step of performing a second 
transformation step comprises: 
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<i> determining an effective range associated with 
each study marker spot member of said set of study marker 
spot members, 

<ii) determining an attraction pairing relationship 
5 between a particular study marker spot member and a 
particular unknown study spot member, said attraction 
pairing relationship being determined utilizing said 
effective range as determined for said particular study 
marker spot member* 

1U (iii) determining positional coordinates in said 

reference coordinate system of said particular unknown study 
spot member by adjusting original scan coordinates by shift 
amounts equivalent to transformation shift amounts of said 
particular study marker spot member that resulted from said 
"first transformation step, and 

<iv> repeating said <ii) and <iii) steps for all 
unknown study spot members. 

8. A method of matching protein spot patterns in a 
plurality of two-dimensional gel electrophoretograms, as 
recited in claim 6, wherein: 

said step of identifying said spot members as marker 
spot members further includes assigning isoelectric focusing 
CPU, and molecular weight <MW> dimensional separation 
values to said marker spot members; and 
25 said method further including, after performing said 

second transformation step, interpolating said assigned PI 
and MW separation values to said marker spot members to 
determine PI and MW values for said unknown study patterns. 

9. A method of matching protein spot patterns in a 
30 plurality of two-dimensional gel electrophoretograms, as 

recited in claim B, wherein: 

said step of comparing includes the step of verifying 
that potentially matching unknown spots are within an 



E0 
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acceptable vector area -formed by marker spot vectors before 
producing said matching results; and 

said producing matching results includes generating a 
matched and unmatched spot data-files -for said unknown study 
5 spot members • 

10- A method of matching protein spot patterns in a 
plurality o-f two-dimensional gel electrophoretograms, as 
recited in claim 9, wherein: 

said producing matching results further includes 

lO generating a cluster datafile having contradictory matching 
information about certain ones of said unknown study spot 
members that match certain other ones of said unknown spot 
members in a contradictory manner; and 

said method further includes resolving said 

15 contradictory matching information, producing spot matching 
results void of said contradictory matching information and 
updating matching results as required to improve accuracy 
and efficiency of the matching task. 

11. A method of matching protein spot patterns in a 
20 plurality of two— d i mens ional gel electrophoretograms, as 
recited in claim lO , wherein: 

said unmatched spot datafile comprises a unique 
datafile having said unknown study spot members with 
coordinates that do not match with any other spot members. 

£5 IS. A method of matching protein spot patterns in a 
plurality of two-dimensional gel electrophoretograms, as 
recited in claim lO, wherein; 

said verification step includes a juxtaposition 
comparison that involves moving a constructed unknown spot's 

30 vector formed from a pair of said potentially matching 
unknown spots towards a pair of marker spot's vectors such 
that their tails have a common point for determining whether 
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said unknown spot's vector is within said acceptable area 
-formed by said pair of marker spot's vectors to verify that 
said potentially matching unknown spots indeed match; and 

said matched spot datafile comprises an exact match 
spot datafile for spot members having been manipulated by 
said verification step. 

13. A method of matching protein spot patterns in a 
plurality of two-dimensional gel electrophoretograms, as 
recited in claim 10, wherein: 

said matched spots datafile comprises a composite 
datafile for producing a pseudo spot pattern representing 
the matched and unmatched spots. 



1**. A method of matching protein spot patterns in a 
plurality of two— dimensional gel electrophoretograms, as 
recited in claim 9* wherein said producing matching results 
further includes: 

generating a single reference gel based matching spot 
dataf i le; 

generating in background a multiple reference gel based 
matching spot datafile having contradictory spot matching 
data and resolving said contradiction and producing a 
composite matching spot datafile; and 

comparing said single reference gel based matching spot 
datafile with said composite matching spot datafile and 
further generating a potential mis— matched spot members 
datafile and a potential matched spot members datafile for 
improving accuracy of said spot matching results. 



15. A method of matching protein spot patterns in a 
plurality of two-dimensional gel electrophoretograms, as 
30 recited in claim 6, wherein 2 

said logistic data and said spot physical 

characteristics data in said scan coordinate system 
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comprises : 

gel record number, spot name , image (gel) name, x 
coordinates, y coordinates, integrated intensity and area. 

16- A method o-f matching protein spot patterns in a 

5 plurality of two-dimensional gel elec trophoretograms (2— D 
gels), said method comprising the steps of: 

< a ) image scanning each of said plurality of 2— D gels 
and producing corresponding spot data files for each 2— D gel 
that contain logistic data and spot physical characteristics 
10 data in a scan coordinate system; 

(b) identifying spot members in each of said plurality 
of 2— D gels that bound investigative spot patterns and 
designating said identified spot members as marker spot 
members; 

15 (c) designating one of said plurality of 2-D gels as a 

reference 2— D gel, said reference 2— D gel having an 
associated reference spot data file including a set of 
reference marker spot members and unknown reference spot 
members having respective coordinates in said scan 

E0 coordinate system; 

(d) designating at least one of remaining ones of said 
plurality of 2-D gels as a study 2-D gel, said study 2-D gel 
having an associated study spot data file including a set of 
study marker spot members and unknown study spot members 

25 having respective coordinates in said scan coordinate 
system; 

(e) performing a first transformation step that 
transforms positional coordinates of said set of reference 
marker spot members, said unknown reference spot members and 

30 each member of said set of study marker spot members from 
said scan coordinate system to a reference coordinate 
system, said first transformation step resulting in each 
member of said set of study marker spot members being in a 
registered relationship with a corresponding member of said 
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set of reference marker spot members; 

<f ) performing a second transformation step that 

transforms positional coordinates of each of said unknown 
study spot members from said scan coordinate system to said 
5 reference coordinate system, said second transformation 
compr isi ng : 

(i) determining an effective range associated 
with each study marker spot member of said set of study 
marker spot members* 

lO (ii) determining an attraction pairing 

relationship between a particular study marker spot member 
and a particular unknown study spot member , said attraction 
pairing relationship being determined utilizing said 
effective range as determined for said particular study 

15 marker spot member T 

<iii) determining positional coordinates in said 
reference coordinate system of said particular unknown study 
spot member by adjusting original scan coordinates by shift 
amounts equivalent to transformation shift amounts of said 

20 particular study marker spot member that resulted from said 
first transformation step, and inputting said determined 
positional coordinates into at least one adjusted data file 
containing new coordinate information for all spot pattern 
members manipulated by said first and second transformation 

25 steps; 

(iv) repeating said <ii) and < i i i ) steps for all 
unknown study spot members; 

(g > repeating said steps ( d ) , (e) and < f ) on any 

remaining one of said plurality of 2— D gels; and 
30 Ch> comparing coordinates of all spot pattern members 

in said at least one adjusted data file and producing 
results identifying matching spot pattern members. 
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17. A method of matching protein spot patterns in 
plurality of two-dimensional gel electrophoretograms , as 
recited in claim 16, wherein: 

said step of comparing includes the step of verifying 
5 that potentially matching unknown spots are within an 
acceptable vector area formed by marker spot vectors before 
producing said matching results; 

said producing matching results includes generating a 
matched and unmatched spot datafiles for said unknown study 

lO spot members? 

said method further includes generating a cluster 
datafile having contradictory matching information about 
certain ones of said unknown study spot members that match 
certain other ones of said unknown spot members in a 
15 contradictory manner? and 

resolving said contradictory matching information and 
producing spot matching results void of said contradictory 
matching information. 

18. A method of matching protein spot patterns in a 
20 plurality of two-dimensional gel electrophoretograms, as 
recited in claim 17, wherein: 

said unmatched spot datafile comprises a unique 
datafile having said unknown study spot members with 
coordinates that do not match with any other spot members. 

S5 19. A method of matching protein spot patterns in a 
plurality of two-dimensional gel electrophoretograms, as 
recited in claim 17, wherein: 

said verification step includes a juxtaposition 
comparison that involves moving a constructed unknown spot's 

30 vector formed from a pair of said potentially matching 
unknown spots towards a pair of marker spot's vectors such 
that their tails have a common point for determining whether 
said unknown spot's vector is within said acceptable area 
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formed by said pair of marker spot's vectors to verify that 
said potentially matching unknown spots indeed match; and 

said matched spot datafile comprises an exact match 
spot datafile for spot members having been manipulated by 
5 said verification step. 

EO- A method of matching protein spot patterns in a 
plurality of two-dimensional gel electrophoretograms, as 
recited in claim 17, wherein: 

said matched spots data-file comprises a composite 
L0 datafile for producing a pseudo spot pattern representing 
matched and unmatched spots - 

El- A method of matching protein spot patterns in a 
plurality of two-dimensional gel electrophoretograms, as 
recited in claim 16, wherein: 
15 said step of comparing includes the step of verifying 

that potentially matching unknown spots are within an 
acceptable vector area formed by marker spot vectors before 
producing said matching results, 

said producing matching results includes generating a 
SO single reference gel based matching spot datafile and 
generating a multiple reference gel based matching spot 
datafile having contradictory spot matching data and 
resolving said contradiction and producing a composite 
matching spot datafile; and 
25 comparing said single reference gel based matching spot 

datafile with said multiple reference gel based datafile and 
further generating a potential mis— matched spot members 
data-file and a potential matched spot members datafile for 
improving accuracy of said spot matching results. 
30 EE. A method of matching two-dimensional <E-D> patterns, 
said method comprising the steps of: 

(a) image scanning a plurality oT E-D patterns and 
producing corresponding data files for each scanned E-D 
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pattern and recording logistic data and pattern physical 
characteristics data in a scan coordinate system; 

(b) identifying pattern members in each of said 
scanned 2— D pattern that bound investigative patterns and 

5 designating said identified pattern members as marker 
members ; 

(c) designating one oT said plurality of 2-D pattern 
as a reference 2-D pattern, said reference 2-D pattern 
having a reference pattern member data file which includes a 

10 set of reference marker members and unknown reference 
pattern members each having respective coordinates in said 
scan coordinate system; 

(d) designating at least one of remaining ones of said 
plurality of 2-D patterns as a study 2-D pattern, said study 

15 2-D pattern having a study pattern data file which includes 
a set of study marker members and unknown study pattern 
members each having respective coordinates in said scan 
coordinate system; 

<e> performing at least one transformation step that 

20 transforms positional coordinates of said plurality of 2— D 
patterns from said scan coordinate system to a reference 
coordinate system for minimizing 2-D pattern preparation 
related distortions; 

(f) repeating said steps (d) and (e) on any remaining 

25 one of said plurality of 2-D pattern; 

(h) comparing physical characteristics of all pattern 
members in said 2— D patterns that were manipulated in 
accordance with said steps <d) through (f ) , said comparing 
including verifying that potentially matching unknown spots 

30 are within an acceptable vector area formed by marker spot 
vectors before producing said matching results; and 

(i) producing matching results. 
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